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Page  1,  Section  13,  Abstract'. 

This  project  extended  understanding  of  vortex  behavior  in  oxide  superconductors  in  two 
important  directions.  First,  we  made  discoveries  regarding  the  morphology  of  the  defects  that 
pin  vortices  in  films  very  effectively  and  regarding  the  nature  of  the  phase  transition  that 
separates  the  low-magnetic-field,  low  temperature  vortex  glass  from  the  high-field,  high- 
temperature  vortex  liquid.  In  both  cases,  our  results  challenge  the  status  quo.  Second,  we 
demonstrated  fabrication  of  YBa2Cu307^  films  with  a  variety  of  dopants  for  Cu,  and  we  showed 
that  only  a  few  percent  of  dopant  has  a  dramatic  effect  on  the  superconducting  properties,  in  a 
way  that  might  be  useful  in  some  superconducting  devices. 


Main  body  of  report,  per  AFOSR  amendment  No.  2  dated  03/27/89. 
a)  Objectives. 

Our  primary  scientific  objective  was  to  understand  the  dynamics  of  superconducting 
vortices  created  in  films  of  the  high-temperature  superconductor,  YBa2Cu307.5,  by  application 
of  a  magnetic  field  perpendicular  to  the  film.  "Understanding"  includes  understanding  of  the 
defect  that  pins  vortices  in  films,  and  understanding  the  nature  of  the  phase  transition  from  a 
vortex  glass  to  a  vortex  liquid.  A  natural  corollary  of  this  project  involved  understanding  how 
vortex  behavior  might  be  influenced  by  the  simple  expedient  of  substituting  a  few  percent  of  the 
Cu  atoms  with  other  atoms  like  Ni  and  Zn. 

Technical  objectives  included  constructing  electronics  appropriate  for  measuring  the 
magnitude  and  phase  of  extremely  small  (nQ-mU)  sample  impedances  at  frequencies  from  dc  to 
20MHz  in  magnetic  fields  up  to  8  Tesla.  Another  objective  was  development  of  sample 
fabrication  facilities  for  making  large  area  films  of  doped  YBa2(Cui.xMj307^,  where  M  =  Ni, 
Zn,  Co,  Al,  etc.  The  final  objective  was  to  find  the  concentration  of  dopants  that  optimized 
vortex  pinning. 


b)  Status  of  research  effort. 

All  necessary  apparati  were  operational  in  the  middle  of  the  second  year,  including  a  10 
Tesla  magnet,  high-speed  electronics,  and  software  for  data  analysis.  Sample  films  of 
YBa2(Cui.xM^)307^  were  made  at  OSU  by  codeposition  and  postannealing;  films  made  by  laser 
ablation  were  obtained  from  the  laboratory  of  Michael  McElfresh  at  Purdue  University.  Data 
were  taken  and  analyzed  by  the  graduate  supported  by  the  grant,  Dimitrios  G.  Xenikos,  through 
the  third  year,  at  which  time  he  finished  his  PhD  thesis  and  left  the  lab. 


Extensive  measurements  of  current-voltage  characteristics  were  made  on  pure 
YBa2Cu307^  and  YBa2(Cuo.98Alo.o2)307-5  films  up  to  20MHz  to  explore  vortex  behavior.  The 
results  of  these  measurements  are  described  below.  As  a  peripheral  project,  films  prepared 
using  fabrication  facilities  supported  in  part  by  this  AFOSR  grant  were  used  to  explore  the  effect 
of  dopants  on  the  infrared  response  of  YBa2Cu307^  films;  hence  publication  2  below  cited 
AFOSR  support. 

The  starting  point  for  understanding  the  behavior  of  vortices  in  YBa2Cu307.5  films  is  to 
understand  the  nature  of  the  defect  in  the  crystal  structure  that  pins  vortices.  Surprisingly, 
despite  the  fundamental  importance  of  this  issue,  very  little  is  known  about  it.  We  have  made 
a  substantial  contribution  by  measuring  the  inductance  of  vortices  by  measuring  the  inductance 
of  a  long  narrow  "bridge"  patterned  into  a  film  as  a  function  of  temperature  and  magnetic  field 
perpendicular  to  the  film.  The  vortex  inductance  is  intimately  connected  to  the  pinning  strength 
-  inductance  grows  as  pinning  weakens.  At  low  magnetic  fields,  the  bridge  inductance  grows 
linearly  with  field,  which  means  that  the  total  vortex  inductance  is  simply  linear  in  the  number 
of  vortices.  At  higher  fields,  the  intervortex  repulsion  weakens  vortex  pinning.  It  turns  out  that 
the  dependence  of  pinning  strength  on  vortex  density  is  extremely  sensitive  to  the  nature  of 
pinning  sites.  We  find  that  pinning  is  weakened  only  gradually  with  increasing  vortex  density, 
consistent  with  pinning  sites  being  lines  or  planes,  but  inconsistent  with  point-like  sites,  which 
is  the  accepted  paradigm.  A  manuscript  describing  this  work  is  in  preparation;  results  appear 
in  the  PhD  thesis  of  Dimitrios  G.  Xenikos. 

An  important  scientific  issue  concerns  whether  the  "melting"  vortices  our  of  their  pinning 
sites,  as  T  is  raised,  is  a  glass-like  softening  or  a  true  melting  phase  transition.  We  have  found 
evidence  which  supports  the  phase  transition  view,  and  evidence  which  has  yet  to  be  reconciled 
with  the  phase  transition  view.  Supporting  evidence  is  that  the  "universal  critical  exponents" 
which  describe  the  power  laws  by  which  various  measurable  quantities  depend  on  T,  or  field  B, 
or  frequency  co,  are  unaffected  when  dopants  are  added  to  the  film.  This  is  nontrivial  since  the 
dopants  have  very  strong  effects  on  various  "nonuniversal"  quantities  including  the  temperature 
at  which  the  transition  occurs.  On  the  other  hand,  by  measuring  the  frequency  dependence  of 
various  quantities,  in  addition  to  the  oft-studied  low-frequency  measurements,  we  discovered  that 
the  two  free  parameters  in  the  theory  were  insufficient  to  describe  the  power  law  behavior  of 
all  of  the  measured  quantities.  Theory  has  yet  to  account  for  this  serious  shortcoming.  The 
former  result  has  been  published;  the  latter  is  in  manuscript  form.  Both  results  appear  in  the 
PhD  thesis  of  Dimitrios  G.  Xenikos. 


c)  Publications. 

1.  "Vortex  correlations  in  the  superconducting  transition  of  YBa2(Cuj_^lJfiy^  films  in 
high  magnetic  fields,"  D.G.  Xenikos,  J.-T.  Kim,  and  T.R.  Lemberger,  Phys.  Rev.  B. 
4S,  7742  (1993). 

2.  "Absence  of  a  gap  in  the  reflectance  of  YBa2(Cuj,J^ifi20  films,"  M.J.  Sumner,  J.-T. 
Kim,  and  T.R.  Lemberger,  Phys.  Rev.  B  12248  (1993). 


4.  "Inconsistent  critical  exponents  for  the  vortex-glass  transition,”  D.G.  Xenikos,  A.E. 
Cunningham,  T.R.  Lemberger,  L.  Hou  and  M.W.  McElfresh,  manuscript  in  preparation; 
to  be  submitted  to  Physical  Review  B. 

5.  "Inductance  of  vortices  in  YBa2CuPy^  films,"  D.G.  Xenikos,  A.E.  Cunningham,  T.R. 
Lemberger,  L.  Hou  and  M.W.  McElfresh,  manuscript  in  preparation;  to  be  submitted 
to  Physical  Review  B. 

6.  "Statics  and  dynamics  of  the  mixed  state  in  pure  and  Al-doped  YBa2Cuj0y^  high 
temperature  superconductors,"  D.G.  Xenikos,  Ph.D.  thesis.  The  Ohio  State  University, 
1994  (unpublished). 
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Consultative  functions: 

Informal  contact  with  Dr.  Timothy  Peterson  at  Wright-Patterson  AFB,  Dayton,  OH. 


f)  Discoveries,  etc. 

The  mobility  of  vortices  can  be  greatly  increased  by  addition  of  a  few  percent  of  A1  or 
Zn  in  place  of  Cu.  This  should  greatly  improve  the  amplification  properties  of  Flux- 
Flow  Transistors  under  development  in  other  laboratories. 
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None. 
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CHAPTER  I 
INTRODUCTION 


A.  On  the  issue  of  true  critical  currents  in  high-Tc  superconductors 

The  early  reports  on  the  statics  and  dynamics  of  the  mixed  state  in  the  high- 
tempeiature  superconductors  (HTSC)  carried  a  great  deal  of  skepticism  over  whether 
true  critical  currents  were  attainable  in  these  new  materials.  It  seemed  that  the  high 
temperatures,  over  which  the  superconductivity  now  persisted,  conspired  with  a  few 
anomalous  properties  of  the  materials  (especially  the  short  coherence  length  and  the 
sample  anisotropy)  to  destroy  their  defining  property  of  zero  linear  resistivity  in  the 
superconducting  state  in  the  presence  of  applied  magnetic  field.  In  fact,  power 
dissipation  in  the  mixed  state  was  an  extensively  studied  phenomenon  in  the  case  of 
conventional  superconductors.  The  new  concern,  however,  was  whether  it  were  now 
possible  to  pin  the  magnetic  flux  quanta  on  the  sample  disorder,  prohibiting  their 
mobility  in  a  dc  current  flow  and  thus  diminishing  the  power  dissipation  in  the  high- 
Tc  samples.  A  huge  effort  was  then  launched  to  understand  these  phenomena,  with 
many  groups  around  the  world  competing  for  innovation  and  applicable  results,  our 
own  group  included. 
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The  issue  of  true  criticsl  currents  was  soon  related  to  that  of  the  nature  of  the 
"irreversibility  line"*’^  Ti„(B).  This  is  a  boundary  in  the  phase  diagram  that  separates 
the  high-temperature  reversible  magnetization  region,  from  the  low-T  irreversible 
magnetization  one.  Initially,  the  similarity  of  the  functional  form  of  the  irreversibility 

line 

Tc-Ti„  a  (1) 

in  HTSC  with  the  Almeida-Thouless  line’  in  spin  glasses  suggested  that  magnetic  flux 
in  the  mixed  state  could  be  considered  as  "frozen"  in  some  random  configuration, 
rather  than  participating  in  a  long  range  magnetic  order  as,  for  example,  in  a 
triangular  array Detailed  magnetization  measurements  though  indicated  that  the 
fluxons  in  the  irreversible  region  were  not  frozen,  but  were  slowly  rearranging 
themselves  to  eliminate  magnetic  field  gradients  in  the  sample  at  logarithmic  relaxation 
rates  ("giant  flux  creep").’  The  fact  that  the  critical  state  could  not  survive  in  HTSC 
crystals  provided  strong  evidence  that  the  disorder  in  the  sample  was  not  effective  to 
pin  the  vortices  in  the  mixed  state. 

Transport  measurements  were  also  interpreted  along  the  same  lines.  The  data 
for  the  temperature  dq)endence  of  the  resistivity  p(T)  were  plotted  to  fit  Arrenius 
curves  p  «  e  and  values  for  the  potential  barrier  U(B,T)  (or  a  distribution  of 
barriers)  were  extracted.*-”  The  suggestion  was  that  thermally  activated  vortices  (or 
bundles  of  vortices)  could  overcome  this  energy  barrier,  consistent  with  an  Anderson- 
Kim  picture”-”,  and  subsequently  dissipate  energy.  Finally  Tinkham”,  by  fitting  the 
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resistivity  data  on  YBajCujOy^  (YBCO)  crystals  of  Inui  et  al.*  using  a  simplified 
phenomenological  analogue  of  the  vortex  motion  over  some  reasonable  potential 
barrier**,  went  all  the  way  to  suggest  that  thermal  vortex  fluctuations  dissipating 
energy  would  dominate  at  liquid  Nj  temperature,  even  in  a  room  temperature 
superconductor.  A  true  critical  state*^’*'  could  not  be  obtained  at  useful  temperatures 
inHTSC. 

The  experiment  that  stirred  the  problem  again  providing  a  new  angle  to  see 
things  was  that  of  R.  Koch  et  This  work  showed  that,  below  the  irreversibility 

temperature  (now  called  glass  temperature  T,),  the  linear  dc  resistivity  was  vanishing 
faster  than  the  Arrenius  behavior,  following  a  power  law  p  «  (T-T,)’^*'**,  with  v(z- 
l)=6.5±0-6.  The  data  were  interpreted  as  evidence  for  a  phase  transition  at  T,, 
separating  a  vortex-liquid  phase,  where  thermal  fluctuations  dominate,  from  a  vortex- 
glass  phase,  with  only  short-range  translational  and  long-range  orientational  order  of 
the  Abrikosov  lattice  surviving.  This  phase  transition  had  been  previously  predicted 
by  M.P.A.  Fisher  on  theoretical  grounds.^* 

From  then  on,  a  number  of  different  theoretical  models,  which  deal  with  the 
statics  and  the  dynamical  characteristics  of  the  mixed  state  and  often  introduce  novel 
phases  with  qualitatively  different  correlations  among  the  vortices,  have  been 
proposed.  The  experimental  evidence  for  the  existence  of  a  phase  transition  at  T,  has 
been  also  building  stronger,  with  current-voltage  (I-V)  curves  in  YBCO  films*’  and 
twinned  crystals^  supporting  the  idea  of  a  phase  boundary  associated  with  the 
"irreversibility  line"^  in  these  materials.  The  details,  though,  of  the  pinning 
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mechanism  that  gives  rise  to  the  critical  behavior  are  still  a  matter  of  controversy. 
Fisher  and  coworkers^'-^  model  vortex  behavior  when  vortex  pinning  is  from 
microscopic  quenched  disorder^  and  they  find  a  second  order  phase  transition 
separating  a  vortex-glass  phase  with  vanishing  linear  resistivity  from  a  vortex-fluid 
phase  with  a  finite  low-voltage  resistivity.  However,  the  nature  of  the  phase  transition 
may  be  quite  different  from  that  of  the  vortex  glass  theory,  if  correlated  disorder  in 
extended  defects  is  responsible  for  the  flux  pinning  in  these  samples.^  Since  both 
models  are  consistent  with  measured  exponents  and  scaling  functions  for  V(I),  one 
must  look  to  non-universal  parameters  to  distinguish  them. 

Our  work  will  move  along  two  directions:  First,  we  will  check  the 
universality  of  the  critical  power  laws  in  a  variety  of  doped  (as  well  as  pure)  samples, 
that  would  support  the  conc^t  of  phase  transition  truly  exists  at  T,.  Second,  we  will 
measure  parameters  that  are  non-universal  in  these  samples,  such  as  the  resistivity  and 
the  pinning  strength,  in  an  effort  to  distinguish  which  of  the  many  types  of  disorder 
is  dominating  in  the  vortex  pinning,  and  relate  to  the  appropriate  theoretical  models. 
Before  we  proceed  to  describe  our  measurements,  we  will  put  forward  in  section  B 
of  the  introduction  the  theoretical  tools  that  we  borrow  from  the  work  in  vortex  glass, 
which  we  will  use  throughout  our  analysis.  In  addition,  in  section  C  of  the 
introduction,  we  will  describe  the  contribution  of  our  work  in  the  field,  and  its 
relationship  with  the  rest  of  the  experimental  and  theoretical  advances  in  the  recent 


literature. 
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B.  Vortex  glass  phenomenology-  the  role  of  the  disorder 

In  this  section  we  will  present  the  basics  of  the  phenomenological  approach  to 
the  glass  transition  which  is  developed  by  Daniel  Fisher,  Mathew  Fisher  and  David 
Huse  (FFH).^  Their  analysis  provides  the  means  to  describe  transport  measurements 
close  to  T,  in  terms  of  "vortex-loop"  fluctuations,  and  will  be  used  extensively  in  our 
work.  We  will  see  that,  although  its  principles  cannot  not  be  confirmed  directly  from 
the  experiment,  the  method  offers  an  insight  on  the  role  of  the  disorder  in  the  sample. 
Here,  we  examine  the  Larkin  length,  the  phenomenology  of  the  vortex  loops,  the 
correlation  length  and  time,  and  finally  the  scaling  characteristics  of  the  transition  at 
T 

As  a  starting  point,  the  vortex  is  assumed  to  be  a  string-like  object  with  certain 
elasticity  eo  carrying  magnetic  flux  4>q,  which  tends  to  be  pinned  in  a  sea  of 
surrounding  microscopic  crystal  imperfections.  Then,  as  argued  by  Larkin  and 
Ovchinnikov^,  the  random  forces  due  to  the  weak  core  pinning  (core  ~  |  in  diameter) 
will  distort  the  flux  lattice  at  large  distances.  Let  us  assume  that  a  transport  current 
J  passes  in  a  direction  perpendicular  to  the  applied  field.  Then,  a  Lorentz  force  acts 
on  the  vortex,  a  segment  L|  of  which  may  become  unpinned  and  move  sideways  a 
distance  Lj.  possibly  smaller  than  L|  (see  Fig.  1).  The  gain  in  the  free  energy 
F(Lj.,L|)  by  this  motion  is  J<^oLj.L|.  The  energy  cost,  given  that  the  vortex  is 
approximated  as  a  string  with  elasticity  eo  per  unit  length,  is  CoLj^VLj.  One  writes 

F(Li,L,)  =  eoLi'/L,  -  -  (if'L,)''’ 


(2) 
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The  third  term  in  the  above  expression  is  the  contribution  from  pinning.  In  the  case 
of  random  point  defects,  as  is  assumed  in  FFH,  the  pinning  energy  <^>  of  a  stiff 
vortex  of  any  length  averages  out  to  zero.  On  the  other  hand,  the  fluctuations  of  the 
pinning  energy  do  not  vanish,  and 


<4(Li)>  =  (3) 

Here,  7  is  a  measure  of  the  pinning  strength,  and  f  the  Ginzburg-Landau  coherence 
length.^*  Then,  minimizing  the  energy  (Eq.  2)  with  respect  to  L|  at  Lj.  and  J=0, 
we  get  the  Larkin  pinning  length  for  the  moving  vortex  segment 


On  length  scales  up  to  the  distortions  of  the  vortex  line  is  smaller  than  the 
coherence  length  (L^  <|),  which  means  that  the  vortex  system  will  look  as  a  lattice. 
However,  at  larger  scales,  "dislocation  loops"  (distortions  of  the  triangular  vortex 
array  where  a  vortex  segment  progresses  one  or  several  lattice  spacings  away,  locally 
rearranging  the  vortex  distribution)  of  length  L|  will  grow,  because  their  contribution 
in  the  free  energy  is  also  €oLx^/L|.^  Thus,  at  lengths  larger  than  the  vortex 
lattice  description  breaks  down  as  far  as  the  translational  order  goes.  On  the  other 
hand,  according  to  FFH  the  orientational  order  of  the  lattice  survives.  This  particular 
vortex  phase  has  been  named  "vortex-glass". 
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The  next  principle  of  the  vortex  glass  phenomenology  concerns  the  energy  of 
the  excitations  of  the  vortex  system.  Making  an  analogy  of  the  vortex  line  in 
quenched  disorder  to  a  "directed  polymer  in  a  random  medium”,  FFH  use  the 
numerical  results  from  the  latter  as  an  indication  that  the  transverse  displacement  Lj. 
of  a  vortex-line  segment  of  length  L|  scales  as 

Lx  -  Lf  (5) 

with  an  exponent  f=0.6,  a  value  larger  than  0.5  (which  corresponds  to  thermal 
fluctuations  and  no  pinning).”  Then  a  typical  vortex  excitation  (see  Fig.  1)  can  be 
thought  as  a  "vortex  loop"  consisting  from  a  vortex  segment  in  the  opposite  direction 
at  the  initial  position  of  the  vortex  segment  that  moved,  connected  to  the  final  position 
of  the  moved  vortex  segment.  The  area  of  the  loop  is  then 

S=  LxL,  -  Li  (6) 

with  jc= 1  +  1/f.  According  to  the  model,  the  vortex  loop  can  involve  many  vortices, 
so  that  Lx  can  exceed  the  intervortex  distance:  A  segment  from  a  vortex  can  move 
to  the  position  occupied  from  another  segment  of  the  nearest  vortex  before  it  move 
to  the  next  nearest  neighbor  vortex,  and  so  on.  It  is  not  clear  if  the  excitations  of  a 
multiple  vortex-line  array  are  preferred  in  the  vortex  system  in  the  HTSC  instead  of 
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Figure  1.  "Vortex  loop"  excitation  (continuous  line).  The  vortex  is 
pinned  on  weak  point  disorder. 
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single  vortex  loop  excitations.  Leaving  all  possibilities  open,  FFH  define  a  "stiffness 
coefficient”  Y  for  the  system  of  vortices  involved  in  the  vortex  loop,  so  that  they  can 
postulate  that  the  elastic  energy  of  the  excitation  is 

YLi  /  L|  ~  Y  (7) 

Then,  balancing  this  elastic  energy  against  the  Lorentz  energy  of  the  vortex  excitation 
J«^oLi.L|  -  JLJ,  FFH  obtain  an  equilibrium  value  of  the  transverse  length  of  the 
vortex  loop  =A  (at  zero  temperature),  in  the  presence  of  transport  current: 

A=(Y/J)''«'f-'>  (8) 

This  is  the  size  to  which  a  fluctuation  must  grow  before  the  current  will  drive  it  off 
to  infinite  radius.  Note  that  the  pinning  information  in  this  analysis  is  all  contained 
in  the  exponent  f,  rather  than  with  an  extra  term  in  the  free  energy  expression  as  in 
the  case  of  Larkin  and  Ovchinnikov  (Eq.  2). 

In  the  next  step,  FFH  estimate  the  electric  field  generated  by  fluctuations  of 
size  A  at  a  temperature  T:  If  there  is  an  energy  barrier  associated  with  such 
fluctuation,  which  relates  to  the  equilibrium  length  A  through  an  (unknown)  exponent 
such  that 


Ba  -AA*, 


(9) 
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where  A  is  some  appropriate  constant  that  indicates  the  barrier  height,  then  the 
resulting  electric  field  generated  by  the  transport  current  in  the  sample  is 

E(J)  -  exp[-BA/kTl  -  exp[-(VJ/]  (10) 

where  /i=^/(2/f-l),  and  Jt-  Y(A/kT),  k=Boltzmann  constant  and  T  the  temperature. 
The  linear  resistivity  in  the  vortex  glass  is  identically  zero  (assuming  of  course  that 
Moreover,  the  expression  (8)  gives  a  nonlinear  voltage  which  becomes 
exponentially  small  for  current  densities  J<  <  Jt-  Thus,  the  theory  of  vortex  glass 
predicts  the  existence  of  true  critical  currents  in  the  mixed  state. 

In  the  vicinity  of  the  glass  transition,  the  theory  resorts  to  another  scaling 
argument,  in  order  to  estimate  the  electric  field.  It  involves  the  "vortex  correlation 
length"  $o>  which  is  defined  as  the  range  of  phase  coherence  in  the  lattice.  This 
means  that  vortex  loop  fluctuations  separated  by  distance  less  than  are  correlated. 
The  vortex  correlation  length  is  expected  to  diverge  at  the  transition  temperature  T, 
as  IT-Tjl",  but  the  exponent  v  is  not  estimated.  In  general,  in  an  anisotropic 
material  the  correlation  length  parallel  and  perpendicular  to  the  field  could  be 
different.  However,  the  vortex  glass  theory  assumes  that  they  diverge  with  the  same 
exponent  i».  Furthermore,  there  is  a  second  quantity  assumed  to  diverge  at  the  phase 
boundary:  the  vortex  correlation  time  -$5,  with  z  the  dynamical  critical  exponent, 
estimated  to  be  in  the  range  z=  4  -7.  Above  T„  this  correlation  time  can  be  thought 
as  the  lifetime  of  vortex  fluctuations  due  to  thermal  excitation. 
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Following  the  above  definitions,  FFH  provide  a  units  argument  to  relate  the 
exponents  z,  p  of  the  theory  to  scaling  expressions  of  the  conductivity  around  T,: 
Since  E  =  -3A/dt  and  the  vector  potential  A  enters  the  Ginzburg-Landau  equation  as 
an  inverse  length,  E  presumably  scales  as  l/(lengthxtime).  In  addition,  J— ^3A, 
where  the  energy  density  /  scales  as  (length)'®  and  D  is  the  dimensionality  of  the 
problem,  suggesting  that  J  scales  as  (length)'"®.  Considering  the  length  f  <3  and  the  time 
scale  ia  of  the  model,  one  concludes  to  a  scaling  expression 

E  -  J  (11) 

E±(x)  are  scaling  functions  for  temperatures  above  (+)  and  below  (-)  the  transition 
temperature  T,,  and  should  be  universal  in  the  materials  where  the  vortex-glass 
transition  occurs.  The  theory  determines  the  scaling  functions  in  certain  limits: 

E.(x)  -  exp(-a/x),  when  x-*0,  (12a) 

£+(*)  ~  E.(x)  -  jajgg  (12b) 

where  a  and  0  are  universal  constants.  Then,  at  T=T,  where  (and  thus 

x-*oo),  a  power  law  is  expected  for  the  current-voltage  characteristic: 


£  ^  j(z+lV(D-l) 


(13) 
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It  is  clear  from  the  above  qualitative  units  analysis  that  a  different  assumption 
for  the  length  and  time  scales  that  are  relative  to  the  problem  would  lead  to  different 
combinations  of  the  critical  exponents  p  and  z  in  the  exponents  of  the  power  laws. 
We  will  address  this  issue  in  our  measurements  in  chapter  III. 

C.  The  role  of  the  disorder  in  the  transition-  theory  and  experiment 

The  nature  of  the  disorder  that  pins  the  flux  lines  in  HTSC  is  still  a 
controversial  issue.  As  we  noted  in  the  previous  section  B,  the  vortex-glass  model 
assumes  that  only  microscopic  quenched  disorder  is  relevant  to  the  problem.  FFH 
maintain  that  their  analysis  could  accommodate  mesoscopic  disorder  (such  as  twin 
boundaries)  in  the  order  of  the  intervortex  spacing  length  scale.  However,  it  is  not 
clear  that  the  exponent  0.5  of  the  theory  (and  thus  all  the  critical  exponents  of  the 
transition)  will  be  unaffected  in  that  case.  Furthermore,  the  assumption  that  a  single 
critical  exponent  p  characterizes  the  divergence  of  both  and  at  T,  could  very 
well  not  hold  for  disorder  of  certain  dimensionality  (2-D  grain  boundaries,  1-D  screw 
dislocations,  etc). 

In  Fig.  2  we  show  some  ways  that  the  flux  lines  can  be  deformed  in  the 
sample.  In  the  left,  microscopic  disorder  (in  the  form  of  perhaps  oxygen  vacancies) 
weakly  bends  the  vortices.  In  the  middle,  we  show  flux  lines  strongly  pinned  in  point 
defects,  when  a  Lorentz  force  is  pushing  them  to  the  right.  Another  alternative  can 
be  pinning  in  columnar  defects  (picture  on  the  right),  where  the  vortices  could  move 
by  the  formation  of  kinks  which  nucleate  by  thermal  fluctuation  or  large  current 
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Figure  2.  Vortex  pinning  in  (from  left  to  right)  weak  point  disorder, 
strong  point  pins,  and  columnar  defects.  A  current  perpendicular  to  the 
page  plane  pushes  vortices  to  the  right.  [Ref.  37] 
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density.  Comparing  the  schematic  figs.  1  and  2,  we  see  that,  in  all  cases  of  defects, 
the  vortex  excitations  could  be  described  as  vortex  loops,  although  the  concept  was 
initially  perceived  for  quenched  disorder.  Starting  from  this  purely  phenomenological 
(unproven)  point,  we  will  keep  the  terminology  of  the  vortex-loop  that  we  described 
in  the  previous  section,  but  we  will  be  ready  to  consider  a  different  relationship 
between  its  transverse  {x  parallel-to-field  components  if  needed. 

There  are  currently  two  theories,  in  addition  to  the  vortex  glass  model,  which 
embody  specific  types  of  defects  and  predict  a  low-temperature  glass  phase  in  the 
mixed  state  of  the  HTSC:  the  collective  pinning  theory  and  the  Bose  glass  theory. 
Contrary  to  the  vortex  glass  model,  where  all  the  information  about  pinning  is 
included  in  the  exponents  f  and  (see  eqs.  5,  10)  through  the  analogue  of  the 
"directed  polymer",  the  collective  pinning  theory  is  a  statistical  solution  of  the 
problem  of  a  collection  of  vortices  in  the  presence  of  point  disorder.  The  theory 
generalizes  the  free  energy  expression  of  Eq.  2,  by  treating  the  flux  lattice  as  an 
elastic  medium,  and  derives  the  exponent  /i  of  the  glass  phase.  It  treats  as  special 
cases  the  statics  and  dynamics  of  a  single  vortex  and  that  of  a  collection  of  vortices 
of  various  sizes  with  respect  to  the  Larkin  length,  and  describes  various  regions  in  the 
H-T  plane  that  may  or  may  not  be  accessible  in  HTSC. 

The  Bose  glass  model,  on  the  other  hand,  considers  vortex  pinning  columnar 
defects,  and  predicts  that  ?|  «  which  leads  to  different  scaling  expressions  in  the 
flux  fluid  and  glassy  phase.  In  particular,  Eq.  11  of  vortex  glass  model  becomes 
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-  £^(J0oe|fx/kT)  (14) 

Then,  one  can  find  a  relationship  between  the  vortex  glass  exponents  p,  z  and  those 
of  the  Bose  glass  p',  z',  by  comparing  with  Eq.  9: 


p’=  1pI7> 

(15a) 

z'=(3z+l)/2 

(15b) 

Furthermore,  the  critical  exponent  in  the  glass  phase  for  the  Bose  glass  is  /*  =  1/3.  We 
will  discuss  the  success  of  the  proposed  theories  to  describe  our  data  in  chapters  III 
and  rV  of  this  work. 

Experiments  testing  the  universality  of  the  critical  scaling  in  a  variety  of 
samples  quickly  followed  the  findings  of  Koch  et  The  direction  of  the  research 
moved  to  either  look  at  samples  of  different  anisotropic  character  (such  as 
Bi2Sr2CaCu20j,  TlSr2Ca2Cu309,  YBCO)  -thus  changing  the  physiognomy  of  the 
vortex,  or  look  at  different  disorder  characteristics  (ex.,  point  or  columnar  defects)  - 
thus  changing  the  size  and  the  strength  of  the  vortex  pinning  sites.  As  the  quality  of 
the  samples  improved,  and  as  the  disorder  becomes  better  characterized,  these 
experiments  revealed  a  wide  variety  of  phenomena  associated  with  T,,  and  led  to 
improvements  of  the  theories. 

Although  there  is  clear  evidence  for  critical  scaling  in  a  variety  of  samples, 
there  is  some  diversity  in  the  published  results.  For  example,  the  self  consistency  of 
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the  critical  scaling  in  YBCO  microtwinned  crystals  was  checked  by  the  AT&T  group” 
and  the  group  at  Caltech”  ”,  but  led  to  different  critical  exponents:  z=5,  p=1.6,  and 
2=s3^  1.2  respectively.  The  scaling  functions  between  crystals  and  films  were  also 

different.  This  was  indirect  proof  that  the  disorder  could  alter  the  nature  of  the 
transition  in  a  qualitative  manner,  and  thus,  special  care  should  be  taken  to 
characterize  the  sample  defects.  Then,  high  quality  untwinned  crystals  showed 
evidence  of  z  first  order  transition  at  the  temperature  associated  with  the  irreversibility 
line.^"”  There  is  evidence  that  this  melting  transition  is  distinct  from  the  vortex  glass 
transition,  which  occurs  only  for  disorder  beyond  a  certain  critical  amount.”  ”  This 
transition  has  been  predicted  from  theoretical  considerations  of  the  effect  of  thermal 
fluctuations  on  the  elastic  moduli  of  the  defect-free  flux  lattice.”"**  Kwok  et  at 
studied  the  effect  of  twin  planes  on  the  transport  measurements  when  the  vortex  is  at 
an  angle  with  respect  to  the  plane  of  the  defects.  Finally,  a  number  of  groups  showed 
a  large  enhancement  of  the  critical  currents  and  a  raise  in  Ti„  in  crystals  with 
columnar  defects.**"*^ 

How  well  does  the  theory  describe  the  diverse  experimental  results  and  the  new 
information  of  specific  types  of  disorder  affecting  the  transition?  First  of  all,  the 
existence  of  a  phase  transition  at  T,  is  still  a  matter  of  controversy,  although  it  is 
consistent  with  the  bulk  of  the  experimental  evidence.  Then,  the  scaling  arguments 
in  the  vicinity  of  T,  or  the  exponent  n  extracted  from  the  data  of  the  glass  phase 
cannot  distinguish  between  the  vortex  glass  model  (in  the  FFH  phenomenological 
picture  or  the  collective  pinning  solution”*^*)  and  the  Bose  glass  model.  In  fact,  the 
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collective  pinning  theory  has  not  explicitly  related  the  pinning  details  of  the  sample 
to  a  phase  boundary  in  the  mixed  state,  although  it  is  consistent  with  a  glass  and  a 
liquid  phase  at  (in  general  unspecified)  parts  of  the  phase  diagram.  The  Bose  glass 
theory  is  definitely  very  successful  in  crystals  with  columnar  defects  (parallel  or 
flayed)  induced  by  heavy  ion  bombardment.  However,  the  predicted  dependence  of 
the  phase  transition  temperature  on  the  angle  of  the  magnetic  field  with  respect  of  the 
defects  was  not  shown  in  films.^  Recent  Monte-Carlo  simulations  based  on  a  mean 
field  approach  by  Stroud  and  coworkers  fit  with  no  adjustable  parameters  the  phase 
boundary  at  YBCO,  in  both  clean  samples  and  ones  with  columnar  defects.^’*® 
However,  some  of  the  results  for  the  pinning  strength  that  we  show  in  chapter  IV 
(especially  the  field  dependence  of  the  Labusch  parameter)  cannot  be  incorporated  in 
the  simulations  at  their  present  form.  Finally,  we  should  mention  that  all  the  above 
theoretical  models  are  seen  with  a  large  degree  of  skepticism  from  a  large  part  of  the 
scientific  community,  which  carefully  looks  for  more  moderate  conventional  ideas  to 
explain  the  experimental  results,  without  the  introduction  of  a  vortex  phase  transition 
at  all. 

In  our  work  we  focused  on  films,  where  little  has  been  done.  We  address 
several  key  issues:  a)  Is  there  a  vortex  phase  transition  in  the  films  at  a  field- 
dependent  temperature  T,(B)?  b)  What  is  the  nature  of  the  low  temperature  phase 
(vortex  glass  or  Bose  glass)  and  what  kind  of  disorder  is  responsible  for  the  critical 
scaling?  c)  How  can  one  tune  the  non-universal  parameters  of  the  transition  without 
affecting  its  universal  characteristics? 
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Thus,  for  the  first  time  in  the  literature  we  varied  materials  parameters  by 
doping  the  YBCO  film.  We  verified  that  the  scaling  exponents  and  functions  of  the 
conductivity  around  T,  are  unchanged,  a  strong  argument  for  the  existence  of  a  phase 
transition  at  T,."-**  Then  our  effort  was  coordinated  on  the  study  of  the  statics  and 
dynamics  of  the  two  vortex  phases  at  both  sides  of  the  transition,  and  find  evidence 
that  correlated  disorder  is  responsible  for  the  transition.  We  propose  several  empirical 
formulas  describing  the  complex  conductivity  of  the  vortex  phase,  and  we  estimate  the 
magnitude  of  the  fundamental  quantities  of  the  transition,  namely,  the  vortex 
correlation  length  and  time.  Finally,  we  suggest  alternative  ways  to  model  the  vortex 
excitations  in  the  light  of  the  vortex-glass  and  Bose-glass  theoretical  models  which  are 


available  at  present. 


CHAPTER  II 


EXPERIMENTAL  TECHNIQUES 

A.  Pure  and  Al-doped  YBaiCujO^  films 

Doping  is  widely  used  for  altering  properties  of  the  normal  and  the 
superconducting  state  in  a  continuous  fashion,  in  order  to  draw  relationships  between 
the  fundamental  quantities  and  the  microscopic  information  at  an  atomic  level.  In  our 
experiment,  A1  was  substituted  for  Cu  in  the  YBCO  lattice,  in  order  to  change  the 
pinning  characteristics  of  YBCO  films.  Yet,  the  role  of  pinning  in  many  of  the 
quantities  that  we  measure,  such  as  the  critical  current  density  or  the  transition  width, 
is  still  quite  a  controversial  issue,  and  forbids  a  quantitative  estimation  of  the  degree 
of  change  in  pinning  that  we  obtain  through  doping.  Consequently,  great  emphasis 
in  this  work  is  given  to  identic  the  nature  of  the  sites  in  the  sample  that  are 
responsible  for  vortex  pinning  phenomena.  It  is  therefore  important  to  discuss  the 
fabrication  technique  and  the  procedures  followed  to  characterize  the  films,  in  order 
to  find  quantities  that  are  sensitive  probes  of  pinning. 

The  films  used  in  this  work  are  made  by  either  codeposition  and  ex  situ 
annealing^,  or  by  laser  ablation.  In  particular,  all  the  doped  films  were  made  by  the 
first  method,  developed  by  Jin-Tae  Kim  in  1991 ,  while  some  pure  laser  ablated  YBCO 
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films  were  used  at  later  stages  of  our  experiments.  The  doped  films  were  prepared 
on  (001)  SiTiOj  substrates  by  codeposition  of  Y,  BaF2,Cu  and  A1  in  the  presence  of 
3  X  10^  Toir  of  O2.  Y  was  evaporated  from  an  electron-beam  source,  while  BaF2,Cu 
and  A1  were  evt^rated  from  heated  filament  sources.  The  evaporation  rates  were 
controlled  independently  by  quartz  crystal  deposition  controllers.  The  films  were 
annealed  ex  situ  in  a  tube  furnace  in  flowing  oxygen.  The  tube  furnace  was  ramped 
quickly  up  to  870  “C  at  40  *C/min  and  then  slowly  up  to  900  “C  at  5*C/min  to  avoid 
thermal  overshooting,  left  at  900  ®C  for  1  hour,  cooled  to  500  ®C  at  5®C/min,  left 
at  500  “C  for  1  hour  to  oxygenate  as  fully  as  possible,  and  finally  cooled  slowly  to 
room  temperature.  While  the  sample  is  at  the  highest  temperature,  900  *C,  the 
oxygen  is  bubbled  through  distilled  water  at  room  temperature.  Annealing  below  900 
°C  or  shorter  than  1  hour  does  not  improve  the  sharpness  of  the  transition,  indicating 
that  it  does  not  raise  the  oxygen  content  of  A1  doped  films.  The  method  was 
successfully  used  for  producing  also  hfi-  and  Zn-doped  films. 

X-ray  diffraction  patterns  for  films  with  less  than  11  %  dopant  concentration 
are  the  same  as  for  highly  oriented  pure  Y6C0  films,  with  no  evidence  for  second 
phases.  Furthermore,  evidence  for  the  high  orientation  of  the  ab  planes  parallel  to  the 
film  plane  is  provided  by  scanning-electron  microscopy  (SEM).  As  an  example,  we 
show  in  Fig.  3  SEM  pictures  of  the  surface  of  a  2%  Al,  1000  A-thick  film.  Clearly, 
there  is  no  evidence  of  misoriented  grains  in  the  surface  topography  which  appear 
only  at  Al  concentrations  in  excess  of  4  %.  The  most  prominent  defects  we  observe 
are  holes  about  1000  A  in  diameter,  which  occur  in  groups  of  1-5  clustered  about  1 
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Figure  3.  Scanning  electron  microscope  pictures  of  a  2  at.  %  Al,  1000 
A-thick  film,  made  by  codeposition  and  ex  situ  post  annealing. 
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ftm  apart.  The  clusters  are  20  apart.  These  holes  may  arise  from  the  collective 
release  of  fluorine  gas  from  the  Aims  during  postannealing.  The  similarity  of  our 
results  of  our  films  made  by  either  codeposition  or  laser  ablation  without  such  defects 
indicates  that  they  did  not  play  any  fundamental  role  -for  vortex  pinning-  in  our 
measurements. 

The  mean  field  transition  was  measured  using  the  apparatus  in  Appendix  B. 
In  Fig.  4  we  compare  the  results  in  A1  with  the  Tc  of  Ni-doped  films  made  by  the 
same  technique.  We  find  that  the  Tc  in  the  Ni-doped  films  dips  about  15  K  over  the 
first  5  %  of  dopant  concentration,  contrary  to  Al-doped  films  where  a  plateau  appears. 
As  will  become  evident  in  chapter  III,  this  attractive  feature  of  the  A1  films  is 
instrumental  when  one  compares  measurements  of  the  transition  width  in  the  presence 
of  magnetic  fields  between  samples  with  different  dopant  concentration. 

Two  of  our  films  were  made  by  the  method  of  pulsed  laser  ablation  by  Lifang 
Hou  and  Mike  McElfresh  in  Purdue  University.  The  films  were  pure  YBCO, 
deposited  on  a  heated  StTiOj  and  a  LaAlOs  substrates.  After  the  deposition,  they 
were  left  in  the  evaporation  system  in  an  O2  atmosphere  for  several  hours  to  anneal, 
as  the  tempCTature  was  gradually  reduced.  As  far  as  the  value  of  Tc,  the  sharpness 
of  the  transition  and  the  normal  state  resistivity,  the  laser  ablated  film  on  SrTiOj  was 
indistinguishable  from  the  codeposited  one.  The  film  on  LaA103  had  a  broad  transition 
with  Tc=82  K  and  was  used  for  setup  calibration  purposes. 

Following  their  characterization  by  ac  susceptibility,  the  films  were  patterned 
into  microbridges  suitable  for  4-terminal  measurements  by  standard  photolithography 
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and  chemical  etching  using  dilute  (0.03  %)  nitric  acid.  The  etching  rate  was  about  400 
A/min.  The  strips  were  typically  15-50  pim  wide  and  100-600  long.  To  reduce 
the  contact  resistance,  approximately  400  nm  of  Ag  was  deposited  on  each  pad,  and 
then  the  samples  are  oxygenated  again  in  a  tube  furnace  at  450-500  ’C  for  1  hour. 
Miniature  coaxial  cables  were  attached  by  pressing  indium  on  the  pads  and  then 
pressing  the  wires  into  the  indium.  The  resulting  contact  resistance  was  typically  in 
the  order  of  mfl  or  less. 

Figure  5  shows  the  ab-plane  resistivity  p(T)  for  four  3000  A  films,  which  is 
measured  by  the  usual  ac  technique,  first  published  in  ref.  56.  The  applied  current 
density  through  the  strip  is  about  30  A/cm^  and  its  frequency  10-100  Kz.  The 
resistivity  p(T)  of  the  doped  films  is  metallic  from  Tc  to  room  temperature,  similar 
to  that  of  the  pure  YBCO  film.  A  transition  to  semiconducting  electron  transport  is 
evident  only  in  the  11  at.  %  Al  film  that  contains  misoiiented  grains,  where  the 
resistivity  exhibits  a  plateau  just  above  Tc- 

After  considering  the  information  from  all  characterization  measurements 
described  above,  we  selected  only  YBajCCunAlJjO^j  films  with  x^0.04  for  our 
measurements  on  the  statics  and  dynamics  of  the  mixed  state.  Our  choice  was 
directed  by  the  need  to  ptftuib  the  pinning  characteristics  in  the  film  without  altering 
them  qualitatively.  Thus,  misoriented  grains  or  non-metallic  normal  state  resistivity 
had  to  be  avoided.  At  low  Al  concentrations,  the  SEM  pictures  of  the  film  surface 
are  indistinguishable  from  those  of  pure  YBCO.  The  T^.  decreases  by  only  ~  1 K  per 
at.  %  dopant  concentration.  The  20%-80%  transition  width  measured  by  either  oc 
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Figure  5.  ab-plane  resistivity  of  3000  A-thick,  Al-dopcd  YBCO  films. 
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susceptibility  or  transport  (at  zero  applied  dc  field)  increases  progressively  from  1-2 
K  at  pure  YBCO  to  4<5  K  for  the  4  at.%  film.  Still,  the  above  do  not  guarantee  that 
the  defects  in  the  sample  (in  the  form  of  point  quenched  disorder,  grain  or  twin 
boundaries,  screw  dislocations  etc.)  are  affected  by  doping.  They  provide  us, 
however,  with  a  variety  of  parameters  controlled  by  doping  that  will  be  finally 
associated  to  the  pinning  characteristics  of  the  flux  lattice. 

B.  Experimental  setup 

The  bulk  of  our  work  consists  of  complex  impedance  measurements  on  the 
high-Tc  films  over  a  wide  frequency  range  10  -  2  x  10^  Hz.  In  general  terms,  we  have 
followed  the  ac  lock-in  technique,  though  the  cabling  and  the  general  configuration 
varies  at  different  frequency  ranges.  Specifically,  the  commercial  lock-in  amplifiers 
detect  signals  typically  up  to  100  KHz.  In  the  radio-frequency  (rf)  range  above  this 
crossover,  the  sample  response  is  heterodyned  into  audio  frequency  signals  that  can 
be  detected  by  the  commercial  lock-in.  The  described  apparatus  has  a  sensitivity  of 
10  nV  and  can  resolve  changes  in  die  impedance  of  a  microstrip  corresponding  to 
a  few  perc^t  of  the  magnitude  of  its  geometrical  inductance,  at  excitation  cuixents 
'-'10'^  A.  Below,  we  describe  the  configuration  and  the  specifications  of  our 
apparatus  in  some  detail. 

The  central  ideas  behind  our  design  can  be  summarized  as  follows: 

1)  Flexibility:  Easy  computer-controlled  frequency  sweep  from  10  Hz  to  20  MHz. 
The  necessary  changes  in  the  cabling  at  the  crossover  f =s  100  KHz  should  only  involve 
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the  room  temperature  electronics. 

2)  Capability  to  measure  the  impedance  of  the  microbridge  which  is  in  the  order  of 
mO  and  switches  over  from  purely  inductive  at  vMy  low  temperatures  to  resistive 
when  approaching  Tc. 

3)  Ability  to  record  both  linear  and  nonlinear  current-voltage  characteristics  in  the 
whole  frequency  range. 

4)  Sensitivity  to  resolve  changes  in  the  order  of  a  fraction  of  the  (mO)  physical 
inductance  of  the  strip,  at  low  current  densities,  where  the  sample  response  is  linear 
to  the  excitation  current. 

5)  Wiring:  The  measurement  should  not  depend  on  the  length  of  the  wires 
connecting  the  sample  to  the  room  temperature  electronics,  so  that  the  film  can  be 
placed  in  the  core  of  a  superconducting  magnet.  Any  background  from  parasitic 
pickup  in  the  cables  should  be  reproducible  on  a  day-to-day  basis  and  consist  only  a 
small  fraction  of  the  measured  signal. 

Figure  6  shows  the  design  of  the  pattern  that  was  used  in  the  wideband 
conducUviiy  measurements.  Dark  colored  is  the  film,  while  is  the  etched  film  area. 
Most  of  the  film  is  used  as  a  ground  plate,  connected  to  the  jacket  of  the  current 
bearing  coaxial  cable,  which  in  turn  is  grounded  to  the  sample  holder  in  the  cryostat, 
1-2  cm  away  from  the  substrate.  We  can  easily  estimate  the  physical  inductance  of 
the  strip  in  such  design:  When  current  Irin(wt)  passes  through  the  strip,  a  magnetic 
field  is  produced  B-[/toI/(2irr)]sin(ut),  with  r>w  the  half  width  of  the  strip.  Then 
the  induced  voltage  between  the  two  voltage  leads  is  V  m  3[^gXud/(2T)]cos(wt), 
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integniting  r  from,  say,  r»w/2=20  /itn  to  2  mm,  with  t  the  length  of  the  strip. 
Then,  the  physical  inductance  of  the  strip  is 

4  -  OJSftot.  (16) 

The  fact  that  the  current  density  peaks  at  the  film  edges  does  not  alter  this  result  by 
more  than  a  few  percent.  For  the  sample  patterns  of  chapter  IV,  f  =600  /tm,  thus 
Lj,»0.5-0.6  nH.  The  signal  we  want  to  detect  will  always  ride  on  this  physical 
inductance. 

The  system  described  here  is  based  on  the  4-terminal  configuration,  which  is 
appropriate  for  wide-range  impedance  measurements  of  a  low-grounded 
superconducting  microbridge^  and  is  the  typical  way  to  measure  low-frequency 
current-voltage  (I-V)  characteristics  of  the  films.  The  components  of  the  complex 
impedance  of  the  load  are  measured  by  a  commercial  lock-in,  whose  bandwidth  is 
extended  by  using  a  heterodyning  technique.  The  variation  of  the  (very  small)  load 
is  overcome  by  careful  termination  of  the  coaxial  cables  that  minimizes  reflections, 
pick-up,  and  the  sensitivity  of  the  measurements  to  the  length  of  the  wires. 

In  Fig.  7,  we  show  two  alternative  ways  that  the  sample  is  connected  to  the 
detection  electronics  for  a  4-terminal  impedance  measurement.  The  one  most 
commonly  used  (Fig.  7b)  has  a  large  current-Umiting  resistor  R= 10^-10*  Q  in  series 
with  the  sample,  and  has  very  poor  response  for  measurements  at  f>  1000  Hz.  jThe 
voltage  drop  across  the  smaller  resistor  R'  =  l-10  0  (R>  >R')  -also  connected  in 
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Figure  7.  Two  configurations  for  4-tenninal  I-V  measurements.  Values 
of  R,  R'  are  given  in  the  text. 
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series  with  the  sample-  indicates  the  current  passing  through  the  bridge].  Instead,  the 
configuration  in  Fig.  7a  has  a  50  0  resistor  in  series  with  the  microbridge,  located  in 
the  cryostat  very  close  to  the  sample.  Thus,  the  coaxial  cable  connecting  the  signal 
generator  at  room  temperature  to  the  sample  (through  a  50  0  symmetric  attenuator 
realized  with  the  resistors  R,  R')  is  properly  terminated  when  the  film  is  in  the 
superconducting  state,  which  reassures  the  absence  of  unwanted  end-reflections  of  the 
current  at  the  sample  up  to  frequencies  in  the  MHz  region.  Moreover,  we  found  that 
the  1  0  resistor  in  series  to  the  sample  for  measuring  the  current  can  be  safely 
eliminated,  simplifying  the  setup  considerably:  the  current  passing  through  the  bridge 
can  be  (to  within  ±3%)  by  dividing  the  (attenuated)  voltage  from  the  signal 

generator  by  the  two-terminal  resistance  of  50  Q.  A  note  of  caution:  the  configuration 
in  Fig.  7a  can  be  used  only  in  the  case  when  the  temperature  was  close  to  Tc  so  that 
the  magnitude  of  the  2-terminal  (complex)  sample  impedance  is  much  smaller  than  50 
0  at  the  frequency  of  the  measurement.  Clearly,  contact  resistance  in  the  order  of  - 1 
mO  is  desirable. 

We  now  emphasize  on  the  connection  of  the  voltage  pads  (in  the  cryostat)  to 
the  instrumentation  preamplifier  (at  room  temperature).  The  sample  is  grounded  at 
one  end  to  the  sample  holder  through  a  short  wire.  A  single  coaxial  cable  (- 1.5  m 
long,  running  parallel  inside  the  probe  to  the  one  that  curries  the  current)  connects  the 
voltage  pads  to  a  floating  50  Q  input  of  the  preamplifier.  [Alternatively,  two  coaxial 
cables  were  connected  to  the  voltage  pads,  with  their  jackets  shorted  at  the  grounded 
end  of  the  film,  for  a  differential  measurement.  This  configuration  was  rather 
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involved  in  that  it  required  good  matching  of  the  two  50  0  inputs  of  the  preamplifier, 
and  gave  no  improved  results  compared  to  the  simpler  cabling.]  In  order  to  properly 
terminate  each  of  these  cables  at  both  ends,  we  could  insert  a  50  0  resistor  between 
its  center  conductor  at  the  cryostat  end  and  the  voltage  pad  of  the  film.  However,  no 
significant  improvement  of  the  background  signal  of  the  setup  is  observed  with  such 
an  addition.  Finally,  the  preamplifier  is  based  on  the  Signetics  592  chip  and  can  be 
used  with  either  single-ended  or  differential  input.  Using  standard  topology  we 
achieved  a  500  Hz  -  22  MHz  (-3dB)  bandwidth  with  30  dB  gain.  We  found  that  one 
amplification  stage  suffices  to  provide  the  necessary  isolation  to  the  sample  and 
dominate  the  mixer  losses,  at  the  high  frequency  measurements. 

For  measurements  at  frequencies  500  Hz  -100  KHz,  the  lock-in  is  connected 
directly  to  the  voltage  leads.  At  10-500  Hz,  the  rf  preamplifier  is  substituted  by  a  low 
noise  commercial  model  (such  as  EG&G  113).  Finally,  at  f  >  100  KHz,  the  voltage 
signal  is  heterodyned  to  be  detected  from  the  lock-in,  as  it  is  shown  in  Fig.  8. 
Following  the  typical  heterodyning  scheme,**  we  mix  both  the  amplified  signal  from 
the  sample  and  the  reference  output  of  the  signal  oscillator  with  a  clean  sinusoidal 
signal  from  a  separate  local  oscillator.  The  frequency  of  the  two  oscillators  differs  by 
Af.  The  output  of  the  IF  port  of  each  mixer  is  fed  to  a  low-pass  filter  (with  a  final 
amplification  stage)  that  provides  the  necessary  50  0  termination  for  the  mixers  and 
attenuates  the  upper  sidebands  (of  frequency  higher  than  Af).  In  most  of  our 
measurements,  a  personal  computer  controls  both  oscillators  keeping  Af=5  KHz.  As 
expected  from  the  broad  bandwidth  of  the  double-balanced  mixers  Synergy  CLK-702B 
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Figure  8.  Heterodyning  scheme  used  for  dynamical  measurements  at 
frequency  range  500  Hz-20  MHz. 
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(0.01-125  MHz  at  the  RF/LO  ports  and  dc-125  MHz  in  the  IF  port),  variation  of  the 
frequency  difference  Af  fiom  1  to  50  KHz  gives  no  detectable  difference  in  our 
measurements.  It  only  sets  a  lower  limit  for  the  frequency  of  the  current  passing 
through  the  sample,  as  heterodyning  a  signal  of  frequency  f  typically  requires 
f>>Af. 

In  Fig.  9,  we  show  the  calibration  of  the  apparatus.  First,  the  bandwidth  of 
our  setup  is  determined  by  bypassing  the  sample  and  applying  a  square  wave  of 
known  amplitude  1  mV  directly  to  the  input  of  the  rf  preamplifier.  Its  output  was 
connected  either  directly  to  our  EG&G  lock-in  for  measurements  at  frequencies 
f<120  KHz  (open  triangles)  or  according  to  the  schematics  of  Fig.  7  for 
measurements  at  f=40  KHz  -20  MHz  (open  squares).  In  the  latter  case,  the  rms 
value  of  the  sinusoidal  signal  at  frequency  Af  measured  by  the  lock-in,  was  converted 
to  its  amplitude  value  by  multiplying  it  by  V2.  The  -3  dB  points  are  500  Hz  and  22 
MHz,  essentially  the  bandwidth  of  the  rf  preamplifier  measured  independently  on  a 
digital  scope.  The  reason  we  calibrated  the  apparatus  with  a  square  signal  and  not  a 
sinusoidal  one  is  that  the  latter  is  -in  general-  deformed  from  a  nonlinear  current- 
voltage  characteristic  of  the  system,  leading  to  a  wrong  estimation  of  the  fundamental 
frequency  when  measured  by  the  lock-in.  On  the  other  hand,  a  square  wave  remains 
square,  even  when  it  passes  through  a  nonlinear  element,  and  the  ratio  of  the 
amplitude  of  the  input  to  that  of  the  output  signal  equals  the  ratio  of  their  components 
at  the  fundamental  frequency,  even  in  the  case  of  a  nonlinear  I-V.  The  impedance 
measurements  that  we  will  discuss  in  chapter  IV  are  taken  in  a  low  current,  so  that 
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Figure  9.  Calibration  of  the  electronic  setup  with:  1  mV  reference 
(top);  inductive  microbridge  (middle);  50  0  resistor  at  the  rf  preamp 
input  (bottom). 
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square  wave  and  sine  wave  signals  give  the  same  voltage  at  the  fundamental 
frequency. 

We  then  measured  the  inductive  response  of  a  600  /tm  x  50  ftm  x  2000  A 
superconducting  strip  at  temperature  below  the  mean  field  transition  at  zero  field. 
This  particular  film  was  made  by  laser  ablation  on  LaAlOj  substrate.  It  had  zero- 
resistance  transition  temperature  at  82  K  and  normal  state  resistivity  of  160  /xO  cm  at 
100  K,  indicating  that  this  sample  was  probably  oxygen  deficient.  In  the 
superconducting  state,  the  sample  exhibited  the  expected  purely  inductive  response: 
at  T=50  K  and  1=  10“*  A  square  wave  current  applied,  the  amplitude  of  the  signal  was 
proportional  to  the  frequency  of  the  current  I  (Fig.  9),  and  the  phase  was  90±3‘’ 
relative  to  I.  The  magnitude  of  the  measured  inductance  is  in  excellent  agreement 
(within  ±  10  %)  with  the  previously  estimated  geometrical  inductance  of  the  strip  (see 
Eq.  12).  The  small  deviation  (±20%)  from  the  linearity  of  the  signal  amplitude  at 
f>3  MHz  is  caused  by  a  background  signal  -lO’’  V,  which  is  very  reproducible 
from  day-to-day  for  the  same  sample  (largely  depending  on  the  details  of  the  cabling 
such  as  the  impedance  characteristics  of  the  contacts  on  the  film)  and  is  independent 
of  the  current  passing  through  the  sample.  The  part  of  this  background  that  is  caused 
by  our  lock-in  technique  can  be  measured  separately  from  the  pick-up  from  the  sample 
cabling,  by  setting  a  50  0  resistor  in  the  input  of  the  rf  preamplifier  in  the  place  of 
the  voltage  cable  from  the  sample.  The  solid  line  in  Fig.  9  indicates  that  this  term 
rises  above  3  x  10^  V  only  for  f  >  3  MHz,  with  a  bandwidth  of  1/3  Hz.  This  noise 
figure  provides  enough  sensitivity  to  resolve  changes  in  the  inductive  response  of  the 
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sample  that  correspond  to  a  few  nanovolts.  It  enabled  us  to  measure  the  kinetic 
inductance  of  YBCO  fUms  in  the  presence  of  a  magnetic  field  over  a  wide  temperature 
anri  frequency  range,  which  typically  represents  a  few  percent  of  the  physical 
inductance.  fThe  ratio  between  the  two  quantities  at  B=0  is  XVdw  *•  1  %,  where  X 
is  the  magnerif:  praetration  depth  -typically  1500  A  at  T=0  K  diverging  at  Tc,  and  d 
and  w  are  the  thickness  and  the  width  of  the  strip,  respectively]. 

We  finish  the  discussion  for  our  experimental  setup  with  some  comments  for 
systems  that  have  been  previously  used  or  are  promising  for  improvement  in  the 
future.  Several  months  after  we  constructed  our  apparatus,  we  learnt  about  a  new  50 
0  lock-in  (EG&G  5202).  It  covers  the  frequency  range  0.1-50  MHz,  thus  it  can 
substitute  the  mixer  portion  of  our  setup,  probably  providing  a  noise  figure 
comparable  to  Fig.  9.  Although  at  the  moment  (March  1993)  this  product  is  rather 
expensive  (3-4  times)  compared  to  our  alternative,  it  is  a  clear  sign  that  the  market 
moves  towards  commercializing  radio-frequency  the  radio-frequency  phase  sensitive 
detection. 

Another  possible  improvement  of  our  setup  may  be  accomplished  with  a  better 
video  ampUfier  in  the  rf  preamplification  stage,  in  place  of  the  Signetics  592.  The 
only  good  choice  available  today  could  be  the  INA103  (InV/v^Hz)  from  Burr-Brown, 
but  with  relatively  similar  bandwidth  characteristics  to  the  chip  we  used. 

Last  but  not  least,  we  should  note  that  measurements  of  vortex  dynamics  in  the 
resistive  transition  at  high  magnetic  fields  are  subject  to  specific  limitations  of  the 
detection  instrument  that  couples  to  the  superconducting  sample.  It  became  apparent 
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in  the  work  of  Olsson  et  a/”,  for  instance,  that,  if  standard  cables  have  to  be  used 
with  the  detection  instrument  (lock-in  amplifier  or  network  analyzer),  the  size  of  the 
cryostat  may  be  severely  restricted,  preventing  the  use  of  superconducting  magnets. 
On  the  other  hand,  permanent  magnets  operating  at  room  temperature  typically  reach 
magnetic  fields  B--0.5  T,  a  value  that  is  in  the  limit  where  dynamical  phenomena 
from  correlated  vortices  can  be  observed.**  Our  setup  avoids  this  difficulty,  since  it 
is  not  sensitive  to  the  length  of  the  semirigid  50  0  cables  connecting  the  film 
(cryostat)  to  the  rf  preamplification  stage  (at  room  temperature).  Our  sample  was 
placed  in  a  vacuum  probe  in  the  core  of  a  9  T  superconducting  magnet,  where  the 
temperature  was  regulated  with  stability  better  than  0. 1  K.  Our  data  of  the  kinetic 
inductance  in  Tesla  fields  below  the  glass  transition  are  among  the  first  in  the 


literature. 
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A.  Universal  characteristics  at  T^T, 

Is  it  possible  to  justify  the  existence  of  a  phase  transition  in  the  vortex  system 
in  superconducting  films  through  only  transport  measurements?  It  is  fair  to  say  that 
unambiguous  proof  for  the  phase  transition  can  not  be  obtained.  There  are  at  least 
two  reasons  for  this:  First,  the  film  geometry  itself  prohibits  test  the  critical  scaling 
predictions  at  temperatures  very  close  to  T,,  where  the  -diverging-  vortex  correlation 
length  exceeds  the  thickness  of  the  sample.  Thus,  conventional  models"  with 
sufficient  adjustable  parameters  may  provide  at  least  semiquantitative  fits  to  the 
transport  data  that  claim  to  provide  evidence  for  the  transition^.  Second,  the  defect 
structure  in  the  films  that  give  rise  to  the  phase  transition  is  extremely  complicated, 
so  that  different  theoretical  models  for  the  transition  (vortex  glass  or  Bose-glass)  could 
be  relevant  to  the  problem.  The  proper  way  to  resolve  this  is  for  one  to  start  with 
samples  with  no  defects,  and  then  introduce  defects  of  a  certain  type  (ex.,  either  point 
defects  by  electron  irradiation  Or  columnar  defects  with  ion  bombardment)  and  study 
the  properties  of  the  transition.  However,  this  can  be  done  only  on  HTSC  single 
crystals.  The  films,  therefore,  provide  us  with  many  challenges:  Since  we  cannot 
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overlook  their  potential  for  technological  applications  we  should  carefully  record  the 
properties  in  their  mixed  state.  We  have  to  consider  all  the  (sometimes 
circumstantial)  evidence  for  the  existence  of  the  phase  transition  and  see  if  they  have 
universal  character  in  many  vortex  pinning  environments.  We  should  finally 
determine  ways  to  control  the  non-universal  quantities  in  the  mixed  state  in  order  to 
determine  the  appropriate  theoretical  model  for  the  transition.  The  way  the  present 
work  addresses  these  issues  is  by  doping  the  YBCO  films  and  by  measuring  the 
current-voltage  characteristics  at  low  frequencies  ("statics")  and  the  film  impedance 
over  a  wide  range  of  frequencies  ("dynamics")  associated  with  the  magnitude  of  the 
vortex  correlation  time  in  the  transition. 

In  this  chapter  we  will  show  evidence  for  a  phase  boundary  in  the  mixed  state, 
drawn  from  current-voltage  characteristics  of  YBa2(Cu,.jAlj307^  films  [0<  x  <0.03], 
in  fields  0.5T  <  B  <  6T.  The  critical  exponents  and  scaling  functions  of  the 
transition  are  independent  of  A1  concentration.  Then,  we  will  determine  the 
characteristic  length  scale  for  vortex  correlations  that  diverges  at  the  transition 
temperature.  We  will  also  show  that  critical  scaling  persists  over  60%  of  the  resistive 
transition  (which  widens  from  3  to  30  K  as  doping  concentration  and  magnetic  field 
increase)  and  breaks  down  when  the  correlation  length  decreases  below  500-1250  A, 
a  value  that  does  not  depend  on  the  intervortex  spacing.  Finally,  we  will  discuss  our 
results  in  the  light  of  the  phenomenological  vortex-loop  picture,  and  suggest  that 
pinning  in  as-made  pure  and  doped  films  of  YBajCoaO?^  is  dominated  by  extended 
rather  than  point  defects. 
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The  films  that  we  report  on  in  this  chapter  are  made  by  co-deposition  of  Y, 
BaFj,  Cu,  and  Al,  with  an  ex  situ  postanneal  in  oxygen.**  They  are  lOOOA  thick, 
single-phase,  and  highly  oriented  with  the  c-axis  perpendicular  to  the  substrate.**  The 
zero-field  transition  temperature  measured  by  ac  susceptibility  in  these  films  is 
suppressed  at  only  1.2  YJ2X.%  of  Al,  in  accordance  with  the  (thicker)  films  of  Fig. 
3,  a  relatively  weak  suppression  compared  to  other  dopants  like  Ni  and  Zn.**  *°  As 
usual,  the  films  are  patterned  for  4-wire  measurements,  with  the  strip  between  the  two 
voltage  leads  now  l(X);tm  x  42ftm.  The  normal  state  resistivity  p(T)  is  metallic,  and 
p(95  K)  increases  from  110  pQcm  for  undoped  YBCO  films  to  150  and  2(X)  pOcm  in 
the  2%  and  3%  Al  doped  films,  respectively. 

Fig.  10  shows  I-V  curves  for  the  2  at.  %  Al  film  at  B=4T  parallel  to  c-axis  at 
temperatures  evenly  spaced  by  approximately  1  K  from  54.0  to  73.4  K  (all  logarithms 
in  this  work  are  at  base  10).  These  I-V’s  are  typical  of  all  films  at  all  fields  studied. 
Following  standard  procedure  in  the  literature,  we  define  the  glass  temperature, 
T,=62±l  K,  to  be  the  highest  temperature  at  which  the  I-V  fits  to  a  single  power 
law,  within  our  voltage  resolution  (50  nV-1  mV).  All  I-Vs  above  T,  exhibit 
systematic  upward  curvature,  thus  restricting  the  uncertainty  of  T,  to  ±1  K.  Fig.  11 
summarizes  T,(B)  [or,  B,(T)]  for  films  with  0%,  2%,  and  3%  Al.  Note  that  T,  for 
the  2%  and  3%  Al  films  is  much  lower  than  for  undoped  YBCO.  Surprisingly, 
doping  at  this  level  weakens  pinning  and  enhances  flux  wandering.  However,  as  we 
will  see  below,  it  does  not  change  the  scaling  functions  and  exponents  of  the 


transition. 
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Figure  10. 1-V  curves  on  a  2  at.  %  A1  film,  at  B=4  Tesla.  The  clashed 
line  indicates  the  transition  at  T,=62  ±  1  K.  Curves  are  taken 
j^roximately  1  K  apart. 
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In  Fig.  12,  we  demonstrate  how  the  I-Vs  above  T,  from  different  samples 
collapse  on  the  same  scaling  curve.  Here,  we  use  the  data  from  the  0%  and  3% 
sample,  at  B=4  T.  The  two  fitting  parameters  are  the  low-voltage  linear  resistivity 
p(T),  and  the  crossover  current  density  Jo’TO  separating  linear  from  nonlinear 
conductivity.  We  arbitrarily  chose  Jo"-(T)  to  be  the  current  density  where  E/(Jp)  =  1.5. 
A  different  choice  for  would  Jo^(T)  would  only  move  the  lines  parallel  to  the  x-axis, 
but  not  affect  the  quality  of  the  fit  which  is  better  than  ±5%. 

The  critical  exponents  of  the  transition  are  also  independent  of  doping.  The 
theory  of  vortex  glass  suggest  ways  to  calculate  z  and  y  from  power  laws  extracted 
from  the  I-Vs.  First,  from  Eq.  13  we  see  that,  V(T=T,)  «  given  that  the 

dimensionality  of  the  problem  is  D=3.  Thus,  the  slope  of  the  I-V  at  T=Tg  (in  a 
logarithmic  plot)  gives  directly  the  critical  exponent  z.  Second,  the  theory  predicts 
p  a  (T-T^'<*‘‘>  and  «  (T-T,)*'.  In  Fig.  13  we  show  typical  resistivity  data  from 
the  3%  film.  Fitting  the  low-temperature  dependence  with  the  predicted  law  we 
extract  the  values  of  v{z-\).  Finally,  the  values  of  extracted  from  the  fitting  of  the 
I-Vs  on  the  universal  scaling  curve  as  in  Fig.  12,  lead  to  the  estimation  of  l-v.  The 
fits  of  V(T=T^,  p(T)  and  Jo‘"(T)  to  the  predicted  power  laws  provide  another  self- 
consistent  way  to  define  T,,  that  leads  to  the  same  value  of  Tj  (within  2  K)  with  its 
operational  definition  from  the  curvature  of  the  I-Vs.  In  Table  1  we  show  the  values 
for  the  exponents  i'(z-l)  and  2p,  as  well  as  (z+l)/2  from  these  fits,  which  are  in 
excellent  agreement  with  the  values  reported  for  YBCO  films*’’^**  and  twinned 
crystals^. 
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Theoretical  models  (other  than  vortex  glass)  would  predict  similar  scaling  laws 
with  different  expressions  of  the  exponents  in  terms  of  p  and  z.  For  example,  we 
mentioned  in  the  introduction  the  suggestions  of  the  Bose*glass  model,  which  is 
proposed  for  samples  with  columnar  defects  rather  than  microscopic  disorder.^  In 
any  case,  the  fact  that  scaling  occurs  with  the  same  exponents  and  scaling  functions, 
even  as  the  width  of  the  resistive  transition  increases  from  3  K  to  30  K  with 
increasing  doping  level  and  field,  supports  the  existence  of  a  phase  transition  at  Tg(B) 
and  our  use  of  appropriate  models  to  interpret  our  data. 

B.  Nonuniversal  characteristics  at 

Despite  the  success  of  the  models  which  embody  the  idea  of  a  vortex  phase 
transition  to  describe  the  scaling  of  the  I-V  curves,  no  consideration  has  been  given 
to  the  interpretation  of  the  nonuniversal  parameters  of  the  transition.  These  include 
the  dependence  of  the  phase  boundary,  the  magnitude  of  the  vortex  correlation  length 
and  time,  and  the  width  of  the  transition  where  scaling  laws  hold,  to  fundamental 
quantities  of  the  supoconducting  state  such  as  the  penetration  and  the  coherence 
length,  or  the  details  of  the  pinning  structure  in  films.  In  this  section  we  address 
these  issues. 

We  first  comment  on  the  reduction  of  Tg  in  the  doped  samples.  For  reference, 
our  2%  film  showed  T,=62±0.5  K  at  B=4T,  to  be  compared  to  Tg=76.3  K  for  our 
undoped  film  and  Tg=79.5  K  for  a  laser  ablated  film  in  ref.  19.  No  theoretical  model 
model  predicts  how  the  glass  transition  depends  on  materials  parameters,  with  the 
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exception  of  the  numerical  simulations  of  Stroud  and  coworkers  for  the  case  of  pure 
YBCO^’*^.  In  general  though,  several  authors  propose  that  T,  may  lie  inside  the 
fluctuation  region  that  is  estimated  from  the  "Ginzburg  criterion". This  criterion 
marks  the  point  at  which  the  system  crosses  over  from  mean-field  to  fluctuation- 
dominated  critical  behavior  as  Tc  is  approached  from  below.  Then,  from  ref  23 
(FFH),  page  138  we  get: 

ATc/Tc  =  (xcA.k/7At)'  (17) 

ATc  is  the  width  of  the  fluctuation-dominated  regime  below  the  superconducting 
transition  Tc=85-90  K,  it=X/foLssl(X)  is  the  Ginzburg-Landau  parameter, 
7  s  (m^/me)*'^  0,2  the  anisotropy  parameter  associated  with  different  effective  masses 
in  the  ab-planes  and  along  the  c-axis,  c,=0.4  a  constant  predicted  by  rough  theoretical 
estimations,  and  a  thermal  length  At=  2  x  10*/Tc  A.  The  zero-temperature 
penetration  depth  is  X„=1600A  and  2800A  in  our  pure  YBCO  and  our  3%  A1  film 
respectively,  measured  by  oc-susceptibility  [E.R.  Ulm  et  dl.  to  be  published]. 
Substituting  in  the  above  equation  we  find  ATc/Tc= 0.15-0.20  for  pure  YBCO,  and 
therefore  T,(B)  lies  inside  the  fluctuation-dominated  regime  for  B<6T.  In  3%  Al, 
ATc  is  3  times  larger,  consistent  with  even  the  widest  transition  we  see  (30K  at 
B=6T).  The  increase  of  Xo  is  common  in  doped  samples  [also  shown  for  Ni-doped 
YBCO:  E.R.  Ulm,  Jin-Tae  Kim  and  T.R.  Lemberger,  to  be  published  in  Physica  B] 
and  0-depleted  YBCO  samples  [W.C.  Lee  and  D.M.  Ginsberg,  Phys.  Rev.  B  44, 
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2815  (1991)]. 

Then,  could  one  predict  the  magnitude  of  the  reduction  of  T,  with  doping  given 
the  increase  in  the  penetration  depth  Xo  due  to  doping?  In  general  terms,  there  should 
be  a  "Lindeman  criterion"  relating  T,  to  parameters  such  as  the  vortex  pinning  energy, 
the  Ginzburg-Landau  coherence  length,  the  penetration  depth,  etc.  In  particular,  the 
«iergy  gain  per  unit  length  from  a  vortex  core  pinned  in  a  defect 

€=  <^o*/(64x^X^  (18) 

decreases  by  a  factor  of  3  with  3  at%  doping,  given  that  the  zero-temperature  Xo,  as 
mentioned  above,  increases  to  2800  A  relative  to  1600 A  of  pure  YBCO.  However, 
T,  decreases  by  only  20%.  It  is  possible  that  A1  introduces  defects  in  the  films  that 
tend  to  raise  T,,  like  screw  dislocations  in  the  crystallites  or  edge  dislocations  in  the 
grain  boundaries.  Such  mechanism  may  dominate  the  reduction  of  e  at  lower  A1 
concentrations  (x  <  0.01)  and  consequently  raise  the  T,  above  its  value  in  pure 
YBCO. 

We  are  now  ready  to  proceed  beyond  scaling  arguments,  to  estimate  the 
magnitude  of  the  vortex  correlation  length  ^o(T)  in  a  model-independent  fashion.  In 
particular,  we  are  interested  in  estimating  $q  at  the  temperature  T  where  the  scaling 
breaks  down,  and  associating  its  value  $*=$o(T*)  to  some  important  mesoscopic 
length  scale  for  vortex  correlations.  We  define  T*  as  the  temperature  where  p(T) 
5%  below  the  value  predicted  by  the  scaling  power  law  p  a  (T-T,)'<*-’».  To 


51 


illustrate.  Fig.  14  shows  p(T)/(T/T,  -  !)'<*•*>  vs.  p(T)  for  the  3%  A1  doped  film  at  B= 
2,  4,  and  6T.  Dividing  out  the  strongest  T-dependence  of  p  emphasizes  the  deviation 
from  the  power  law.  The  value  of  T*  and  thus  our  conclusions  are  not  particularly 
sensitive  to  our  criterion  of  5%,  which  determines  T*  to  ±1  K.  Interestingly,  the 
critical  scaling  persists  over  60±5%  of  the  width  of  the  resistive  transition,  i.e. ,  (T*- 
Tg)/(Tc-T^*»0.60,  even  though  the  width  varies  from  3  K  to  30  K.  Here,  we  take 
the  onset  of  the  resistive  transition  T®  to  be  where  p(T)  in  the  applied  magnetic  field 
is  about  75%  of  the  extrapolated  normal  state  resistivity.  This  corresponds  to  a 
suppression  in  the  mean-field  transition  temperature  of  —0.6  K/T,  about  the  same  as 
in  pure  YBCO.  The  significance  of  this  "universal"  behavior  for  T*  is  unknown. 
From  the  applied  point  of  view  though,  doped  samples  may  thus  be  proved  a  useful 
tool  in  case  a  large  critical  temperature  regime  is  required,  as  their  broader  resistive 
transitions  correspond  to  wider  critical  regions. 

We  now  estimate  ^oCT)  using  the  crossover  current  density  which  we 
define  from  the  scaling  forms  as  the  current  density  J  where  E/pJ=1.5  (see  Fig.  12). 
Since  all  I-Vs  in  the  critical  regime  collapse  on  the  same  scaling  curve,  other  choices 
of  this  ratio  between  1  and  3  do  not  affect  our  results.  We  use  the  conventional 

expression^’* 


-  kT 


(19) 


that  does  not  depend  on  the  particular  model  for  the  phase  transition,  to  determine 
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lo"  »  (20) 

Physically,  this  relation  means  that  nonlinear  behavior  occurs  when  the  current  density 
is  large  enough  to  produce  translations  of  correlated  vortex  segments  of  length  over 
distance  |j.  perpendicular  to  the  current,  as  readily  as  thermal  fluctuations  do.  At 
translations  of  correlated  segments  are  thermally  induced,  while  at 
Lorentz-force-induced  translations  dominate.  Since  oe  (T-T,)^',  the  size  of  vortex 
fluctuations  diverges  at  T,  as  consistent  with  an  increasing  number  of 

vortices  becoming  correlated  and  moving  coherently  as  T-*T,.  To  simplify  the 
notation,  we  define  r  =  (T/T,  -  1),  r*  «  (T*/Tg  -  1),  and  write 


fo=  e*(rVr)' 


(21) 


Fig.  15  shows  $a‘^=  Jo-^CD^o/kT  vs.  t/t*  for  a  3%  A1  film  for  B=  2  and  6T.  The 
fact  that  the  two  data  sets  overlap  means  that  both  the  magnitude  and  scaling  law  for 
|o(T)  are  independent  of  B,  and  hence  independent  of  the  average  intervortex  spacing. 
The  intercept  at  r/T’=l  for  this  film  yields  =  5 15  ±40  A,  which  exceeds  somewhat 
the  intervortex  spacing  ay  =  (<^o/B)'^  —  190A  and  320A  for  B=6T  and  2T, 
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respectively.  If  were  proportional  to  ay,  then  the  2  T  data  (filled  symbols  in  Fig. 
15)  would  fall  on  the  dashed  line.  Taking  into  account  the  uncertainties  of  our 
estimation  in  T,,  T*  and  the  exponent  i'(z-l),  we  find 

r  a  B®*®,  (23) 

i.e. ,  an  extremely  weak  dependence  of  on  the  applied  field.  This  is  the  main  result 
of  this  work  and  it  holds  for  all  the  samples  we  studied,  and  for  the  laser  ablated  films 
of  ref.  19  (see  Fig.  16).  Table  1  records  the  values  of  obtained  for  each  film.  A1 
doping  up  to  3  at.  %  has  only  a  small  15%  effect  on  the  magnitude  of  |  .  On  the 
other  hand,  evidently  film  morphology  affects  |*.  In  a  25(X)  A-thick  YBCO  film  with 
a  sharp  transition,  we  found  f  *  =  1240±60  A.  In  an  unusually  resistive  2500A-thick 
YBCO  film  with  a  wider  superconducting  transition  but  the  same  scaling  exponents 
we  found  |’=  730±40A.  Analyzing  the  measured  V(T)  for  B=  2,  3  and  4T  from 
Koch  et  al.^^  on  laser-ablated  4006A-thick  YBCO  films  (Fig.  16),  we  find  ^  = 
1100  ±200  A,  independent  of  the  magnetic  field. 

What  is  the  significance  of  this  result?  First  we  consider  the  case  that 
strongly  depends  on  the  specifics  of  the  behavior  above  the  transition  7* .  This  is 
suggested  by  M.  Salamon  et  al*^  who  note  that  T*  may  mark  the  crossover  between 
two  regimes  with  different  scaling  characteristics.  They  propose  that,  in  the 
temperature  region  centered  on  T^  and  estimated  by  the  Ginzburg  criterion,  a  3-D 
X-Y  critical  scaling  model"  describes  the  film  conductivity.  This  argument  is 


Figure  17.  Scaled  conductivity  a(T,H)  data  according  to  X-Y  model  in 
ref.  61  from  a  pure  and  a  3  at.  %  Al  film  do  not  collapse  onto  a 
universal  line.  t=l-T/Tc. 
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consistent  with  our  observation  that  both  boundaries  T*(B)  and  T,(B)  scale  with  field 
with  the  same  power  exponent  (which  in  principle  can  be  related  to  a  3-D  critical 
exponent 


(Tc-T*)  «  (Tc-T^  oe  v^-0.1  (24) 

However,  an  attempt  to  scale  our  p(T)  data  in  the  vicinity  of  T  to  the  expected  form*' 

0(7, H)  H'"=  (25) 

[H»B,  o=  1/p  is  the  low-current  ohmic  conductivity,  f  «  (Tc-T)"'  is  the  coherence 
length]  leads  to  non-universal  scaling  functions  jKCOH*'^,  strongly  dependent  on 
dopant  concentration.  This  is  demonstrated  in  Fig.  17,  where  the  left-hand-side  of  Eq. 
25  is  plotted  against  =(1-T/Tc)/H®  ‘'*’  «  [f(T)H''^-‘"'.  Even  if  we  consider 

an  increase  of  the  magnitude  of  {  to  accommodate  the  broadening  of  the  transition 
with  doping,  still  the  two  curves  would  move  relative  to  each  other  parallel  to  the  x- 
axis  of  the  plot,  and  they  would  not  fall  on  the  same  scaling  curve  as  Eq.  25  suggests. 
During  my  presentation  in  LT-20  conference.  Dr.  A.  Koshelev  mentioned  that  the 
feilure  of  this  particular  model  by  Salamon  et  d  could  mean  that  the  problem  is 
considerably  more  complicated  that  is  previously  thought,  rather  than  that  the  3-D  XY 
scaling  model  is  completely  inapplicable.  Definitely,  more  theoretical  work  has  to  be 
done  to  clarify  the  issue. 
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We  now  approach  the  problem  from  the  phenomenological  point  of  view  of  the 
vortex  loops,  and  consider  the  magnitude  of  in  terms  of  vortex  correlations.  In  any 
model  for  the  phase  transition  at  T,,  the  critical  response  must  involve  correlated 
motion  of  multiple  vortices.®  Hence,  the  correlation  length  fo  (or  at  least  f  x)  should 
be  larger  than  the  average  intervortex  spacing  in  the  critical  region.  Models  in  which 
pinning  is  hrom  quenched  disorder  with  an  average  spacing  between  pinning  centers 
much  smaller  than  the  intervortex  spacing  lead  to  the  conclusion  that  scaling  should 
break  down  when  the  correlation  length  is  roughly  the  intervortex  spacing  (a  field 
dependent  quantity).  In  the  vortex  glass  theory  in  particular,  the  spatial  extent  of 
phase  coherence  of  the  superconducting  order  parameter  is  assumed  to  exceed  the 
range  of  translational  order  of  the  flux  lattice.  Thus,  the  critical  scaling  is  expected 
to  break  down  when  decreases  below  the  Larkin-Ovchinnikov  length^,  which  is  at 
most  a  few  fiux-lattice  spacings.”  Again,  should  be  quite  sensitive  to  B  if  a  bundle 
of  vortices  are  involved  in  the  problem,  and  sensitive  to  T,  as  shown  from  Eq.  4  even 
in  the  case  of  a  single  vortex,  in  sharp  contrast  to  our  observations. 

A  natural  interpretation  of  our  result  is  that  reflects  the  underlying  pinning 
structure  in  the  films.  Quite  generally,  models  for  critical  scaling  in  samples  where 
correlated  disorder  dominates  pinning,  require  the  transverse  displacement  Sx  of 
liberated  vortex  segments  to  be  at  least  the  distance  between  neighboring  macroscopic 
defects  in  the  sample^,  and  thus  lead  to  a  field  independent  value  for  at  the 
temperature  where  critical  scaling  breaks  down.  Equally  probable  sites  for  flux 
pinning  in  our  films  are  grain  boundaries,  twin  planes  perpendicular  to  the  substrate, 
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and  columnar  defects  along  the  c-crystallographic  axis.  A  typical  grain  size  of  500 
A  in  YBCO  films  has  been  measured  by  TEM®,  a  value  close  to  our  f  .  Pinning  in 
twin  planes  is  also  shown  to  dominate  the  regime  above  the  T,  in  twinned  YBCO 
crystals.”’^  Finally,  assuming  that  only  columnar  defects  such  as  screw  dislocations 
in  crystallites  or  edge  dislocations  in  grain  boundaries  are  responsible  for  the  pinning, 
a  Bose-glass"  calculation  leads  to  a  reasonable  ds=  200  A  for  the  distance  between 
adjacent  defects."  Although  our  experiment  cannot  decide  for  the  nature  of  the 
extended  defects  responsible  for  the  critical  scaling,  it  suggests  that  their  structure  is 
dictated  by  the  morphology  of  the  substrate  and/or  misfits  in  the  substrate-film 
interface  rather  than  the  dopant  concentration,  since  A1  doping  has  only  a  15%  effect 
on  I*.  Furthermore,  the  variation  of  with  film  thickness  in  undoped  films  suggests 
that  sample  defects  change  somewhat  when  the  film  thickness  exceeds  the  typical  grain 
size.  It  remains  to  be  understood  how  A1  doping  weakens  the  pinning  properties  of 

the  extended  defects. 

Below,  we  give  the  argument  of  the  Bose-glass  theory  leading  to  the  estimation 
of  the  distance  between  correlated  defects  mentioned  in  the  previous  paragraph,  that 
was  communicated  to  us  by  David  Nelson.  In  ref.  24,  Nelson  and  Vinokur  show  in 
detail  how  the  correlation  length  relates  to  the  distance  d  between  columnar 

defects: 


{,=d(T/TiO' 


(26) 
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with  Tm*  60  K  a  characteristic  temperature  in  the  model,  f  |  iVD,,  where  1 1  and 
are  the  parallel  and  perpendicular  to  the  field  correlation  lengths.  Do  =*3x10^  cm 
is  a  diffusion  length,  suggestive  of  how  far  out  from  the  applied  field  direction  does 
a  vortex  segment  deviate  relative  to  its  length  f|.  In  our  2500A  film  we  have 
(f  |{i)‘'*=1250A  at  T=T*,  and  thus  d=200  A.  In  our  1000  A  -thick  films,  d-250 
A. 

As  a  final  point  in  the  chapter,  we  discuss  the  striking  result  that  the  values  of 
p**p(T=T*)  are  clustered  around  22  pO  cm  for  all  our  samples  and  at  all  fields  (inset 
in  Fig.  15).  Laser  ablated  films  also  show  p**25  pO  cm  at  fields  B=2-4T.‘  A 
sample-independent  p*  sets  a  strong  constraint  for  the  magnitude  of  the  vortex 
correlation  length  ^  «  (T-T,)  '  and  time  t  «  with  v,  z  the  static  and  dynamical 
critical  exponent  of  the  transition.  In  the  framework  of  the  vortex  glass  model,  p  oc 
(T-t^kx-i)  anji  thus  p  oe  f /t.  Therefore,  at  T=T*,  the  ratio  |/t  is  independent  (similar 
to  p*)  on  the  applied  magnetic  field,  the  width  of  the  transition  or  the  normal  state 
resistivity  of  the  sample.  The  magnitude  of  $  and  t  can  be  measured  independently 
by  and  ac-”  conductivity  measurements  respectively,  and  provide  quantitative 
checks  for  the  theoretical  models  that  describe  the  critical  response  of  the  vortex 
system. 

Although  above  we  used  a  vortex-glass  argument,  the  fact  that  the  B- 
independent  values  of  p*,  set  a  strong  constraint  for  the  magnitude  of  the  relaxation 
time  t=  t*($o/|*)*,  associated  with  the  critical  slowing  down,  is  fairly  robust.  We 
estimate  t*  by  considering  a  vortex  excitation  of  area  formed  when  current  density 
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passes  through  the  sample  at  temperature  T* .  The  energy  gained  in  the  system  by 
the  generation  of  the  excitation  is  Jo''"^o{**“kT .  As  this  excitation  relaxes,  it 
energy  at  a  rate  J  J*e  (fx,  with  J  and  e  being  the  local  current  density  and 
local  electric  field  caused  by  the  excitation,  integrated  over  volume.  Assuming  that 
the  integral  is  significant  only  at  the  core  of  the  vortex  of  radius  Jol,  and  that 

c*p-V(ayioL)  (27) 

is  the  macroscopic  electric  field  E»p*Jo*  increased  by  the  geometric  factor  we 

obtain 

f  J-edfx  -  (28) 

Letting  t*  be  the  time  for  the  fluctuation  to  dissipate  its  energy  of  kT,  we  find 

f-CrOfavloL^-p*)  (29) 

Olsson  et  find  from  the  complex  impedance  of  a  2500A-thick  YBCO  film  at 
B=0.55  T  that  critical  scaling  persists  for  t  ^  2xl0^s.  Using  values  for  ^  and  p* 
from  our  2500  A  pure  YBCO  film,  and  ^ql^IOO  A  at  T  =88  K,  we  estimate 
t*»3xl0^  s,  in  reasonable  agreement  with  the  experiment.  As  B  increases  and  T 
moves  away  from  Tc(B),  inflated  values  of  |gl  are  avoided  and  the  system  should 
become  slower  with  t*  *  a,'*.  One  should  note  however,  that  a  different 
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phenomenological  picture  for  the  vortex  excitations  could  lead  to  a  new  combination 
of  the  relevant  length  scales  (including  possibly  an  anisotropic  fo)  ***  estimation 
of  t*.  The  issue  may  be  resolved  by  dynamical  measurements  in  an  extended  region 
of  magnetic  fields  and  film  thicknesses,  and  will  be  partially  addressed  in  the 
following  chapter.  These  measurements  have  the  potential  to  estimate  directly  the 
relaxation  time  of  the  vortex  excitations,  which  refiects  the  dynamics  of  the  fiux  lines 
as  they  cross  and  recombine  in  the  sample. 


TABLE  liCritical  scaling  parameters  of  YBa2(CU|.,Alj307^  films.  The  uncertainties 
represent  possible  systematic  errors  in  the  calculations. 


CHAPTER  IV 


DYNAMICS 


A.  Glass  phase  at  T<  <T, 

In  the  previous  chapter  we  examined  the  scaling  characteristics  in  the  vortex 
liquid  state  above  the  transition  and  identified  the  parameters  that  are  universal  and 
those  which  vary  with  doping.  We  now  proceed  to  study  the  pinning  characteristics 
of  the  glass  vortex  phase  below  T,.  Our  goal  is  to  determine  experimentally  the 
strength  of  the  pinning  force  exerted  on  the  vortices  due  to  the  sample  disorder,  so 
that  a  detailed  description  of  the  interaction  between  magnetic  flux  and  defects  can  be 
achieved.  Currently,  there  are  two  theoretical  models  that  start  from  specific 
considerations  about  the  sample  disorder  and  lead  to  a  melting  transition:  the 
collective  pinning  theory”’*',  which  treats  the  flux  lattice  as  an  elastic  medium  in  weak 
quenched  disorder  and  concludes  that  a  vortex  glass  phase  exists  at  low  temperatures, 
and  the  Bose-glass  model”  for  vortices  strongly  pinned  on  correlated  disorder,  which 
leads  to  a  Bose-glass  vortex  phase  below  T,.  Furthermore,  recent  numerical 
simulations  that  incorporate  correlated  disorder  are  also  promising  in  calculating  the 
vortex  pinning  strength  and  are  consistent  with  flux  melting  at  The 

physiognomy  and  the  strength  of  the  vortex  pinning  sites  in  the  sample  assumed  in 


64 


65 


these  models  are  quite  diverse  and  should  be  determined  experimentally,  so  that  a 
correct  description  for  the  observed  phenomena  be  reached. 

Direct  observation  of  pinning  effects  can  be  obtained  from  the  electrodynamic 
response  of  the  samples,  Le.,  the  study  of  the  complex  response  functions  such  as 
impedance  and  susceptibility.**^  In  particular,  the  nc-conductivity  of  the  HTSC  in 
the  presence  of  dc  magnetic  field  B  |c  has  shown  direct  evidence  for  critical  slowing 
down  above  the  melting  transition”*^,  and  evidence  for  "intrinsic  pinning",  Le., 
pinning  of  a  vortices  on  the  ab-plane  when  the  magnetic  field  is  applied  at  small 
angles  with  the  ab-crystallographic  plane**.  In  our  work,  we  expand  the  technique  to 
study  the  film  inductive  response  at  the  temperatures  below  the  onset  of  dc  dissipation 
[at  T=T,(B)  estimated  independently  from  current-voltage  (I-V)  curves].  For  the  first 
time  we  find  a  phenomenological  relationship  between  the  pinning  strength  of  the 
vortices  (those  being  macroscopic  objects)  in  the  glass  phase  and  the  inductive 
response  of  the  supercurrents  at  zero  field  (which  directly  relates  to  the  magnetic 
penetration  depth,  a  fundamental  superconducting  parameter).  Thus,  we  find  that:  1) 
Each  vortex  in  the  sample  feels  the  same  pinning  potential,  2)  the  effective  restoring 
force  constant  due  to  pinning  is  inversely  proportional  to  the  square  of  the 
temperature-dependent  of  the  magnetic  penetration  depth,  and  3)  the  sample 
admittance  drops  to  zero  over  a  few  degrees  Kelvin  in  the  vicinity  of  T,,  indicating 
that  the  sample  disorder  becomes  ineffective  to  pin  the  flux  lattice  above  that 
temperature.  In  addition,  we  find  the  ac-impedance  of  the  sample  in  the  vicinity  of 
the  melting  temp^ature  to  be  consistent  with  critical  slowing  down  of  the  vortex 
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system.  We  discuss  our  results  in  the  light  of  the  different  types  of  disorder  that  are 
present  in  the  films. 

The  epitaxial  YBCO  thin  film  studied  in  this  part  of  our  work  was  deposited 
by  pulsed  laser  ablation  onto  a  heated  SiTiO,  substrate.  The  film  was  d=2500±400 
A  thick,  with  its  c-axis  perpendicular  to  the  film  plane,  and  it  exhibited  a  sharp 
transition  (the  20% -80%  transition  width  was  less  than  2  K).  The  zero-temperature 
penetration  depth  of  the  ab-planes  was  X(0)=1730±100  A  measured  with  the  2-coil 
method,  in  good  agreement  with  that  of  high  quality  crystals.**  Following  our 
established  procedure,  the  film  was  patterned  into  a  4-contact  microbridge  of  width 
w=35  urn  and  length  f=600  /xm  by  usual  photolithographic  techniques  and  wet 
etching  in  dilute  nitric  acid.  The  contacts  were  defined  by  ion  milling  (to  clean  the 
film  surface)  and  subsequent  Ag  deposition  (0.5  fim  thick).  The  zero-resistance 
temperature  of  the  film  at  B=0  was  88.5  K. 

Measurements  were  taken  in  the  transport  microbridge  using  a  square  pulse 
current  of  amplitude  I  (much  less  than  the  critical  current  of  the  strip)  and  frequency 
f=10^  -  2x10*  Hz  and  a  lock-in  to  measure  the  amplitude  of  the  signal  V  from  the 
sample.  The  current  was  kept  small  compared  to  the  critical  current  of  the  strip,  so 
that  the  measured  oc-voltage  V  was  linear  in  I  at  all  temperatures  of  interest  (elastic 
limit**).  For  frequencies  f=50  KHz-  20  MHz,  it  was  necessary  to  heterodyne  the 
sample  response  into  audio  frequency  in  order  to  be  detected  by  the  lock-in,  using  the 
apparatus  in  chapter  II.  Special  care  has  been  taken  to  eliminate  reflections  in  the 
transport  lines  from  load  mismatch  by  inserting  50  0  resistors  in  series  with  the 
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sample  (which  is  essentially  a  short,  as  its  complex  impedance  never  exceeds  10  mfl 
in  our  measurements).  In  that  way,  we  managed  to  reduce  parasitic  signals  to  a 
magnitude  2-10  times  smaller  than  that  from  the  physical  inductance  of  the  sample, 
throughout  the  frequency  region  this  inductance  was  measurable,  and  at  excitation 
currents  as  low  as  ~  30  /xA. 

In  Fig.  18  we  plot  the  magnitude  jZ^I  (inset:  phase  of  the  impedance  Z,, 
=  (V/I)dw/f  =  1  Zb  I  e**  measured  accross  the  microbridge  corrected  by  the  geometrical 
factor  dw/f  so  that  we  can  relate  Z^  to  microscopic  sample  parameters,  as  a  function 
of  magnetic  field  B  applied  parallel  to  the  c-axis,  at  f  =  1.  MHz.  The  average  current 
density  I/(dw)  used  in  the  measurement  was  6x10®  Am'^,  /.e.,  as  low  as  1-3  orders 
of  magnitude  smaller  than  the  critical  current.  [The  actual  current  density  is  expected 
to  peak  at  the  film  edges,  and  this  will  be  taken  into  consideration  below.]  From  the 
figure  we  see  that  the  bridge  impedance  Zj,  can  be  separated  into  three  components 

Zb  =  Zp  +  Zk  +  Z,  (30) 

with  distinct  physical  meaning.  First,  a  term 

Zp=a)Lp(dw/f)  (31) 

due  to  the  magnetic  flux  coupled  into  the  voltage  circuit  (dashed  line),  can  be 
estimated  to  about  +20%  accuracy  from  its  geometrical  dimensions  (Eq.  17)  as 
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fjLf^  •  0.5-0.6  nH.  This  term  is  treated  as  an  adjustable  parameter  in  our  analysis, 
and  it  is  important  to  consider  its  T-  and  B-dependence.  Such  dependences  can  arise 
from  a  change  of  the  distribution  of  the  current  in  the  strip,  which  may  alter  the 
magnetic  flux  coupled  between  the  voltage  leads  of  the  microbridge.  This  current 
distribution  depends  on  the  T-dependence  magnetic  penetration  depth  X(T)  and  the 
local  vortex  density  in  the  film.  However,  X(T)  changes  by  at  most  a  factor  of  2  in 
the  temperature  region  below  80K  where  most  of  our  data  are  taken,  and  the  vortex 
density  is  reasonably  uniform  with  an  exception  probably  a  distance  comparable  to  the 
effective  penetration  depth  2\Vd  from  the  film  edges.  A  simple  estimation  of  the  flux 
coupled  to  the  voltage  leads  shows  that  does  not  change  across  the  T-  and  B-range 
of  our  experiments  by  more  than  a  few  percent. 

Given  the  above  considerations  regarding  Lj,,  the  y-intercept  of  the  solid  lines 
in  fig.  18  represents  the  temperature  dependent  zero-field  kinetic  inductance 


Zk=a,MoX'(T)  (32) 

from  the  inertia  of  the  superconducting  carriers  (a}=2Tf  the  angular  frequency,  and 
X(T)  the  magnetic  penetration  depth).  Finally,  the  data  at  B>0  indicate  the  effect  of 
the  driving  current  on  the  vortices.  Driven  by  the  Lorentz  force  JxB,  the  flux  quanta 
oscillate  and  give  rise  to  a  -generally-  complex  conductivity  0y=ayi+ia^2—  1/Zy.  The 
data  in  Fig.  18  and  its  inset  set  a  constraint  for  the  real  and  imaginary  part  of  the 
conductivity:  in  the  region  where  1 2^,  |  a  1/B,  the  sample  response  is  purely  inductive 
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(^=90°  within  our  experimental  error  ±  0.4“).  Conveniently  expressing 

1/[wmoXc2(T,B)]  (33) 

in  terms  of  the  phenomenological  Campbell  penetration  depth  Xc(T,B)  in  a  notation 
similar  to  the  form  of  the  zero-field  kinetic  inductance,  the  above  constraint  suggests 
that,  at  temperatures  where  the  signal  is  purely  inductive  Xc^(T,B)cx  B.  This  is  the 
first  important  result  of  our  dynamics  measurements. 

We  now  study  the  T-dependence  of  Xc(T,B)  by  comparing  it  with  X(T).  The 
data  in  Fig.  19  are  taken  at  temperatures  where  the  sample  response  is  strictly 
inductive.  Using  a  value  of  Zp  which  is  within  10%  of  our  estimated  value,  the  B=0 
data  lead  to  a  magnetic  penetration  depth  X=X(0)/(l-T/Tc)‘^,  where  X(0)= 1730  A  lies 
within  the  uncertainty  of  the  value  measured  independently  by  ac  -susceptibility  in  the 
same  film  before  patterning.  A  change  of  our  assumption  for  Zp  by  ±  20%  does  not 
affect  the  inferred  T-dependence  of  X(T)  «  (1-T/Tc)‘*^  at  T>60  K,  and  adjusts  the 
x-intercept  of  the  regressions  in  Fig.  19  (corresponding  to  Tc)  by  less  than  ±  1  K. 
At  B=1  T,  the  T-dependence  of  the  signal  remains  the  same,  thus 

(X^+Xc^)-'  oc  (Tc-T)  (34) 

Actually,  because  X'^(T)  has  some  downward  curvature  above  60  K,  we  expect  the 
intercept  of  X’^(T)  to  hit  above  the  measured  Tc.  Combining  the  information  from  Eq. 
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34  and  our  previous  result  Xc^ocB,  we  write 

Xc"  =  (B/Bo)  X"  (35) 

with  Bo=0.48±0.05  T.  This  means  that  the  inductance  of  the  vortices  equals  the 
kinetic  inductance  of  the  superconducting  electrons  at  B=0.5  T.  In  general,  we  find 
that  a  single  value  of  Bo  describes  all  our  data  (as  long  as  a^j  =0)  for  0.5  <  B  <  6.5 
T,  and  up  to  our  highest  frequency  f=20  MHz.  Furthermore,  Bq  does  not  change 
in  the  presence  of  a  dc  current  passing  through  the  bridge,  provided  that  it  is  less  than 
the  critical  current.  Finally,  although  we  looked  for  it,  we  could  not  see  any 
hysteresis  in  our  field  data  whether  B  was  increasing  or  decreasing,  or  its  polarity  was 
reversed.  Using  a  time  constant  of  11  s  in  the  lock-in,  possible  systematic  deviations 
from  Eq.  (35)  are  within  our  experimental  random  error.  Eq.  (35)  relates 
quantitatively  the  intrinsic  magnetic  penetration  depth  X(T)  with  a  macroscopic 
quantity  Xc(T)  and  is  one  of  the  central  results  of  this  work. 

We  analyze  the  T-  and  B-dependence  of  our  data  for  T<  <T,,  in  terms  of  the 
Labusch  parameter  which  is  the  curvature  of  the  pinning  potential  at  its 
minimum.’^’*  Driven  by  the  Lorentz  force,  F=^oJsin«t,  each  flux  quantum 
oscillates  over  distance  x=doSinwt  giving  rise  to  an  inductive  term  in  the  measured 
complex  impedance  (rv2'‘=«Bd</»s«t.  Taking  the  phenomenological  standpoint  that 
each  vortex  oscillates  under  a  restoring  force  per  unit  length  ]  F 1  =aLX,  and  using  ^,2'* 


=ti)/ioXc’,  we  find  using  Eq.  35 
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V(T)=B<^o/(MoaL). 

(36) 

(37) 

Since,  at  any  given  temperature,  ai,  is  observed  to  be  constant  for  transport 
currents  up  the  critical  current,  the  elastic  limit  the  width  of  the  potential 

where  our  approximation  for  its  quadradic  form  breaks  down)  is  determined  by  the 
critical  current  density  Jc(T): 

=  <^oJc(T)/aL(T)  (38) 


Since  both  JcCD  and  0^(7)  are  approximately  proportional  to  Tc'T  in  the  YBCO 
films,  we  find  that  f  .i  is  in  the  order  of  10  A  almost  independent  of  temperature,  i.e. , 
it  is  approximately  the  size  of  the  vortex  core  in  the  sample.  Given  the  form  of  Ol, 
the  height  of  the  pinning  potential  at  {*1  is  in  the  order  of  10  meV,  a  value  much 
smaller  than  estimations  for  the  potential  bamer  height  derived  from  magnetic 
relaxation  experiments  in  YBCO  (typically  200  -300  meV)* 

Another  important  consequence  of  Eq.  31  is  that  it  sets  a  semiquantitative  limit 
in  the  range  over  which  the  strength  of  the  pinning  overcomes  the  intervortex 
interactions  and  thus  destroys  the  translational  order  of  the  Abrikosov  lattice  in  the 
YBCO  film.  We  notice  that 
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(39) 

i.e.,  that  the  magnitude  of  the  pinning  restoring  force  constant  a,,  in  the  range  of 
fields  of  our  experiment  is  of  the  same  order  as 

c«=BV(4T/ioX^  (40) 

the  mean-field  shear  modulus  of  the  triangular  lattice^^.  This  result  indicates  that  the 
strength  of  pinning  is  comparable  to  the  energy  needed  to  shear  the  flux  lattice,  and 
is  consistent  with  estimations  from  dc  I-V  curves  at  T<T,  that  the  Larkin  length  in 
YBCO  is  only  a  few  intervortex  spacings.’’  Furthermore,  the  fact  that  the  mismatch 
between  the  Labusch  parameter  and  the  shear  modulus  is  temperature  independent 
(both  are  proportional  to  1/X^)  suggests  that  the  range  over  which  the  flux  lattice  is 
ordered  does  not  change  with  temperature  at  T<T,,  exactly  as  expected  in  a  "frozen" 
vortex  state. 

By  far,  however,  the  most  surprising  result  in  this  work  is  that  field 
independent.  This  underlines  that  each  vortex  is  subjected  to  the  same  restoring 
pinning  force,  rather  than  a  distribution  of  pinning  potentials.  This  fact  excludes  the 
possibility  that  vortices  are  pinned  on  isolated  pinning  sites,  hierarchically  filling  the 
deepest  potential  wells  first  and  then  (as  more  vortices  enter  the  sample  by  the 
increase  of  B)  occupying  weaker  defect  sites.  It  rather  suggests  that  the  number  of 
pinning  sites  in  the  film  greatly  exceeds  the  number  of  vortices.  Then,  in  the  case 
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that  pinning  in  correlated  disorder  dominates,  one  concludes  that  only  macroscopic 
defects  (as  columnar  defects)  of  certain  size  provide  effective  pinning.  According  to 
the  work  of  SdSik  and  Stroud,  the  pinning  strength  of  such  sites  resembles  that  of 
columnar  defects  of  diameter  50  A,  where  the  superconducting  order  parameter  is 
reduced  by  10%.  On  the  other  hand,  pinning  due  to  microscopic  disorder  in  the 
sample  cannot  be  ruled  out  from  our  data.  In  this  eventuality,  each  vortex  core  may 
rest  on  a  large  number  of  defects,  thus  averaging  out  their  strength  and  leading  to  a 
Labusch  parameter  independent  of  the  vortex  density  or  the  position  of  vortices  in  the 
film.  In  this  case,  the  collective  pinning  theory  should  in  principle  describe  our 
experimental  form  of  Ol,  the  magnitude  of  which  could  be  estimated  considering 
either  vortex-core  pinning  or  magnetic  pinning  of  the  fluxons  on  the  point  defects. 


B.  Critical  region  at  T»T, 

We  now  proceed  to  study  the  linear  response  of  the  vortices  at  temperatures 
and  magnetic  fields  where  the  conductivity  of  the  strip  has  a  component  in-phase  with 
the  excitation  current  (<rvi?^0).  In  this  region,  the  sample  conductivity  is  frequency 
dependent,  as  is  demonstrated  in  Fig.  20.  There  we  show  the  reactance,  l/ay2(f,T), 
measured  at  an  applied  field  B=6  T,  after  we  subtracted  the  data  taken  at  the  same 
temperature  but  zero  field.  The  straight  lines  represent  fits  0^,2'*  «  f“,  with  the 
exponent  a  determined  without  adjustable  parameters  from  the  slope  in  the  logarithmic 
plot,  and  shown  separately  in  Fig.  21a.  Such  power  law  behavior  has  been  shown  to 
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f  (Hz) 


Figure  20.  Logarithmic  plot  of  av2**  vs.  frequency,  at  temperatures  66.9,  74.6,  75.3, 
75.9,  and  77.0  K.  The  straight  lines  correspond  to 
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be  consistent  with  the  existence  of  a  phase  transition  of  the  vortex  system.” 

A  reasonable  question  that  arises  is  whether  the  observed  scaling  is 
quantitatively  consistent  with  previous  results.  To  check  it,  we  took  I-V  curves  at 
f=100  Hz  (dc  limit)  and  estimated  the  transition  temperature  T,(B=  6T)  =  73.7±0.7 
K  as  the  lowest  temperature  with  non-zero  linear  resistivity  (using  the  analysis  of 
Refs.  20,  22).  Then,  at  T=T,,  Fig.  21a  shows 

a =0.82  ±0.04,  (41) 

corresponding  to  vortex-glass  dynamical  critical  exponent  z=6±l.  The  same  value 
of  or  is  found  from  analysis  of 


a^,-‘(T=T,)  oc 

(42a) 

^  a  ojosstoo*  (at  T=T,). 

(42b) 

In  addition.  Fig.  22  shows  that  4>(T=T,)=72  ±  4®,  a  value  which  is  predicted 
directly  from  a  Kramers-Kronig  transformation  of  the  power  law  in  Eq.  (29b).™  This 
value  is  in  excellent  agreement  with  previous  work  on  complex  impedance 
measurements  around  T,,  and  is  consistent  with  critical  scaling  of  the  dc 
measurements  (vortex  glass  dynamical  critical  exponent  z=6±l).  In  Fig.  21b  we  go 
beyond  scaling  arguments  to  show  a  novel  characteristic  suggestive  of  a  phase 
transition  at  T,:  a  drop  in  the  admittance  <«)(rv2(T)  =  which  becomes  sharper 
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as  frequency  is  reduced.  At  T<T,,  the  sample  inductance  is  nonzero  (it  follows  Eq. 
(27))  and  proportional  to  f  (so  that  «  «)»  the  admittance  is  a  nonzero  f- 
independent  quantity.  On  the  contrary,  at  T>T,,  Ovz'*  «  w®  with  o<  1.  Therefore, 
the  admittance  drops  to  zero  over  a  T-region  that  decreases  as 

Although  Fig.  21b  is  suggestive  of  such  trend,  finite  size  effects  do  not  allow 
the  transition  to  become  very  narrow:  analysis  of  dc  I-V  curves  (as  in  chapter  HI) 
indicate  that  the  vortex  correlation  length  |g(T)  at  B=6  T  exceeds  the  film  thickness 
in  a  region  of  3-4  K  around  the  glass  temperature  T,.  Thus,  the  transition  cannot  be 
narrower  than  what  is  indicated  in  our  Fig.  21b,  even  if  we  had  enough  resolution  to 

take  inductance  data  at  lower  frequencies. 

What  is  the  significance  of  this  result?  A  drop  in  the  film  admittance  at  zero 
field  in  conventional  superconductors  has  been  correlated  with  a  Kosterlitz-Thouless 
phase  transition.’’’’^  In  the  presence  of  a  net  magnetic  field,  the  film  admittance  is 
related  to  the  perpendicular  to  the  field  component  of  the  helicity  modulus  tensor  by 
the  form  yj.  «  wa^z,  taken  in  the  limit  of  zero  frequency.^  Thus,  the  drop  in 
from  its  frequency-independent  finite  value  (imTvi  “  1/(Mo\:^)  at  T<Tg  to  zero  over  a 
few  degrees  above  T,,  indicates  a  sharp  decrease  of  the  superfluid  density.  This 
effect  is  recently  predicted  by  §45ik  and  Stroud”  using  Monte  Carlo  simulations  on 
disordered  YBCO  samples,  and  is  a  strong  evidence  of  melting  of  the  flux  lattice  at 

T 

1,. 

Summarizing  the  work  in  this  chapter,  we  have  measured  the  field  (0.5  <  B 
^  6.5  T),  temperature  (50  <  T  <  90  K)  and  frequency  (50  KHz  <  f  <  20  MHz) 


Figure  22.  (a)  I-V  curves  plotted  as  log(V/I)  vs.  log(I);  f-dependence  of  the  amplitude 
of  V  (b)  and  the  phase  of  V/I  (c),  with  1=  40  fiA.  Temperatures  are  approximately 

0.6  K  apart. 
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dependence  of  the  complex  conductivity  of  z  YBCO  microstrip.  We  found  a  drop  in 
the  admittance  at  T»Tg  (Tg  was  estimated  independently  from  dc  I- Vs).  Scaling  at 
Tg  is  consistent  with  the  existence  of  a  phase  transition  of  the  vortex  system.  At 
T<Tg  we  extracted  the  restoring  force  that  pins  the  flux  in  the  sample.  Finally,  we 
estimated  the  elastic  limit  where  a  quadratic  pinning  potential  well  describes  our 
dynamics  data,  and  related  it  to  the  possible  types  of  defects  in  the  sample  that  may 
dominate  the  low-temperature  elastic  response  of  the  flux  lattice. 


CHAPTER  V 


CONCLUSIONS 


In  this  work  we  attempted  a  thorough  study  of  vortex  pinning  effects  as  they 
are  realized  in  statics  and  dynamical  response  of  the  flux  lattice  on  transport  currents. 
We  proceeded  in  two  directions:  First,  we  tested  the  universality  of  phenomena  that 
were  previously  observed  in  pure  YBCO  films,  and  are  supposed  to  hold  in  a  variety 
of  high-Tc  materials  with  similar  defects  in  their  crystal  lattice.  To  this  end,  we  used 
films  with  small  concentration  of  A1  substituting  the  Cu  in  the  atomic  cell.  Second, 
we  derived  for  the  first  time  empirical  laws  that  describe  quantitatively  the 
temperature  dependence  of  the  complex  conductivity  in  the  mixed  state,  related 
quantities  in  the  superconducting  transition  as  the  transition  width  and  the  magnitude 
of  the  linear  resistivity,  and  estimated  the  fundamental  quantities  such  as  the  vortex 
correlation  length  and  time. 

In  some  of  the  issues,  our  data  are  strongly  suggestive.  For  example,  the 
question  "are  there  true  critical  currents  in  high  temperature  superconductors?"  is 
answered  with  a  clear  "yes":  All  our  data  are  consistent  with  a  phase  transition  at  a 
temperature  T,(B)  from  a  vortex  phase  with  linear  dc  resistivity  to  a  glass  phase  where 
the  flux  lattice  is  effectively  pinned  in  the  sample  disorder.  In  particular,  we  find  that 
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the  static  nonlinear  I-V  characteristics  obey  power  laws  with  critical  exponents  and 

functions  that  are  independent  of  doping  (and  thus  independent  of  numerous 
quantities  that  characterize  the  films  as  transition  width,  normal  state  resistivity, 
London  penetration  depth  etc.).  For  the  first  time  we  showed  a  drop  in  the  sample 
admittance  at  the  T,(B)  that  is  defined  from  the  I-V  curves.  Finally,  the  critical 
slowing  down  was  directly  observed  in  the  frequency  response  of  the  ac  conductivity 
at  rf  frequencies. 

Parallel  to  clarifying  the  above  issues,  our  data  revealed  a  deep  weakness  of 
theoretical  models  to  relate  the  measured  quantities  (resistivity,  value  of  T,(B),  etc.) 
to  fundamental  parameters  in  the  superconducting  state  as  the  penetration  depth,  or 
to  information  on  the  microscopies  of  the  disorder  in  the  films.  It  is  certainly  quite 
true  that  an  experimental  work,  such  as  the  present  one,  is  always  allowed  to  stop  at 
a  phenomenological  level  of  understanding  of  the  phenomena  at  hand,  attempting  only 
to  describe  them  in  a  self-consistent  manner.  However,  such  fact  can  hardly  let  the 
mind  ease,  at  the  prospect  of  a  better  arrangement  of  ideas  that  we  are  so  certain  that 
will  come.  In  the  meanwhile,  what  I  will  only  add  at  the  end  of  this  manuscript  is  a 
table  of  the  empirical  formulas  resulted  through  this  work,  as  a  monument  -I  hope- 
of  my  extreme  dedication  to  the  experimental  detail  that  conducted  my  way. 
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Table  2.  Summary  of  the  phenomenological  laws  in  this  work. 


1 

Empirical  formulas 

Prefactors 

1  linear  resistivity 

P=P*{_£I5z.)7-o*o.3 

T*-Tg 

p’=22±5  /tOcm 

correlation  length* 

T*-T 

f*=500  -1250A 

correlation  time 

T*-T 

t=t*(  '‘y)7.0±0.3 

T-Tg 

t*-30ns 

Labusch 

parameter 

- 

Bo=0.48±0.05  T 

1  transition  width** 

r-T,-0.6  (Tc“-T,) 

- 

"The  value  of  depends  on  the  doping  level  and  the  film  thickness,  but  is  field 
indq)endent. 

•*  Tc“  denotes  the  mean  field  transition  temperature  when  field  B  is  applied. 


APPENDIX  A 


NONLINEAR  MAGNETIZATION  OF  Y-B-C-0  CRYSTALS 

The  early  experimental  work  -late  nineteen  eighties-  in  YBCO  run  on  two  main 
themes:  First,  the  need  to  record  the  novel  character  of  these  samples,  especially  at 
temperatures  close  to  the  mean  field  transition  where  the  irreversibility  line  was 
identified.  The  result  was  a  steady  stream  of  publications  from  numerous  laboratories 
that  dealt  with  the  peculiarities  of  the  resistive  transition,  in  addition  to  the  question 
of  the  existence  of  true  critical  currents  in  the  samples.  Second,  the  effort  to  produce 
continuously  improved  samples  in  the  form  of  single-phase  powder,  large 
monocrystals  or  thin  films.  Our  lab  quickly  engaged  in  the  effort  in  producing  high- 
quality  single  crystals  and  films. 

Following  the  successful  fabrication  of  crystals  with  the  desired  sharp 
transitions,  we  developed  a  method  to  measure  their  magnetization  in  the  presence  of 
superimposed  ac  and  dc  magnetic  fields  in  the  order  of  a  few  Gauss.  In  particular, 
we  focused  on  the  harmonic  content  of  the  magnetization,  which,  in  YBCO  ceramics, 
has  been  shown  to  be  surprisingly  rich  and  had  raised  questions  about  the  nature  of 
the  critical  state  in  the  high-Tc.  To  our  surprise,  we  realized  that  a  semiquantitative 
understanding  of  the  phenomenon  could  be  obtained  by  considering  a  toy-model  for 
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the  screening  currents  flowing  in  a  resistive  medium,  where  the  resistance  could  be 
borrowed  from  available  results  from  transport  measurements.  Our  dynamics  model 
was  based  on  the  observed  "flux  creep"  resistance,  in  which  the  temperature  T  and 
the  magnetic  field  B  appear  only  in  the  ratio  [The  reader  should  refer 

for  a  detailed  discussion  of  our  results  to  ref.  79,  which  for  technical  reasons  cannot 
be  included  in  this  manuscript.]  Clearly,  more  complicated  theories  should  have  -at 
the  very  least-  the  same  success  in  fitting  the  data,  something  that  was  not 
accomplished,  as  far  as  we  know.  It  would  be  quite  interesting  for  one  to  repeat  these 
measurement  in  crystals  with  well  characterized  disorder  (twin  planes),  which  are 
available  today,  to  confirm  some  of  these  early  findings. 


APPENDIX  B 


AC  SUSCEPTIBILITY  APPARATUS  FOR  MEASURING  THE  TRANSITION 
TEMPERATURE  OF  HIGH-Tc  CRYSTALS,  SINTERED  SAMPLES,  AND 

FILMS 

The  ac  susceptibility  technique  is  a  reliable  technique  that  can  probe  the  onset 
of  the  diamagnetic  response  in  a  superconducting  sample.  As  a  non-contact  technique, 
it  is  non-intrusive,  allowing  the  sample  under  consideration  to  be  used  in  subsequent 
measurements.  Since  it  does  not  require  a  contingent  superconducting  path  between 
two  points  where  contacts  are  attached  to  the  sample  (as  in  a  transport  measurement), 
it  can  measure  the  superconducting  volume  of  the  sample  and  reveal  the  intergrain  or 
intragrain  origin  of  its  diamagnetic  response.  Finally,  by  recording  the  transition 
width,  it  leads  directly  to  information  on  the  London  penetration  depth  (when  mGauss 
excitation  magnetic  fields  are  used)  or  the  properties  of  the  critical  state. 

The  ac  susceptibility  method  can  be  utilized  by  using  3  coils  (one  primary  and 
two  detection  coils  as  in  Appendix  A),  2  coils  (one  excitation  and  one  detection,  as 
it  is  currently  used  in  our  lab  to  measure  the  penetration  depth  in  films),  or  finally  1 
coil  in  a  tank  circuit,  as  is  described  in  this  appendix.*®  Experimental  details  about 
the  various  methods  of  ac  susceptibility  can  be  found  in  ref.  81.  Our  setup,  in 
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particular,  serves  as  an  extremely  flexible  device:  in  the  past  five  years,  it  was  used 
hundreds  of  times  to  characterize  sub-millimeter-size  crystals,  thin  films  of  area  1 
cm^,  and  sintered  samples  of  various  size  used  as  targets  in  laser  ablation  and 


sputtering  of  YBCO  films. 
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We  determine  the  characteristic  length  scale  for  vortex  correlations  in  YBa2(Cui_;5Al;j  )307_5  films 
[05;c  £0.03]  by  studying  their  current-voltage  characteristics  in  fields  0,5  T£R  £6  T.  We  find  Jtiat 
critical  scaling  persists  over  60%  of  the  resistive  transition  (which  widens  from  3  to  30  K  ^  doping  con¬ 
centration  and  field  increase)  and  breaks  down  when  the  correlation  length  is  5(X)— 1250  A,  a  value  that 
does  not  depend  on  the  intervortex  spacing.  This  result  suggests  that  pinmng  in  as-made  films  of 
YBa2Cu307^5  is  dominated  by  extended  rather  than  microscopic  defects. 


Vortex  pinning  in  disorder  high- temperature  supercon¬ 
ductors  has  generated  a  number  of  different  theoretical 
models,  which  deal  with  the  static  and  dynamical  charac¬ 
teristics  of  the  mixed  state  and  often  introduce  phases 
with  qualitatively  different  correlations  among  the  vor¬ 
tices.  There  is  strong  experimental  evidence  from 
current- voltage  (/-K)  curves  in  YBa2Cu307_5  (YBCO) 
films  ^  and  twinned  crystals^  supporting  the  idea  of  a 
phase  boundary  associated  with  the  “irreversibility  line”^ 
in  these  materials.  The  details,  though,  of  the  pining 
mechanism  that  gives  rise  to  the  critical  behavior  are  still 
a  matter  of  controversy.  Fisher  and  co-workers^*^  model 
vortex  behavior  when  vortex  pinning  is  from  microscopic 
quenched  disorder,  and  they  find  a  second-order  phase 
transition  separating  a  vortex-glass  phase  with  vanishing 
linear  resistivity  form  a  vortex-fluid  phase  with  a  finite 
low-voltage  resistivity.  However,  the  nature  of  the  phase 
transition  may  be  quite  different  from  that  of  the  vortex- 
glass  theory,  if  correlated  disorder  in  extended  defects  is 
responsible  for  the  flux  pinning  in  these  samples.  ^  Since 
both  models  are  consistent  with  measured  exponents  and 
scaling  functions  for  F(/),  one  must  look  to  nonuniversal 
parameters  to  distinguish  them. 

In  this  paper,  we  investigate  the  type  of  disorder  that 
pins  vortices  in  the  glass  state  in  high-T^  films  by  examin¬ 
ing  the  length  associated  with  vortex  correlations 
above  the  glass  transition.  We  report  that  critical  scaling 
persists  over  temperatures  where  the  spatial  extent  for 
phase  correlations  of  the  superconducting  order  parame¬ 
ter  exceeds  a  length  (500-1250  A)  that  does  not  de¬ 
pend  on  the  intervortex  separation.  Such  mesoscopic 
length  scale  is  absent  from  models  that  assume  pinning 
only  in  pointlike  disorder,  including  the  vortex-glass 
model.  To  investigate  the  generality  of  this  result,  we 
dope  our  YBCO  films  by  substituting  up  to  3%  of  Cu 
with  Al.  This  substitution  leaves  the  critical  exponents 
and  scaling  functions,  theoretically  universal  quantities, 
unchanged,  although  it  enhances  thermal  fluctuations 
and  reduces  T^{B).  The  length  is  hardly  affected  by 
doping,  and  it  is  independent  of  the  intervortex  separa¬ 
tion. 

The  films  are  made  by  co-deposition  of  Y,  BaF2,  Cu, 


and  Al,  with  an  ex  situ  postanneal  in  oxygen.^  They  are 
1000- A  thick,  single-phase,  and  highly  oriented  with  the  c 
axis  perpendicular  to  the  substrate.^  Doping  suppresses 
the  zero-field  transition  temperature  measured  by  ac  sus¬ 
ceptibility  at  only  1.2  K/at.%  of  Al,  a  relatively  weak 
suppression  compared  to  other  dopants  like  Ni  and  Zn.^’* 
The  films  are  patterned  with  standard  photolithographic 
techniques  for  four-wire  measurements,  with  the  strip  be¬ 
tween  the  two  voltage  leads  typically  100  X  42  /xm.  The 
normal  state  resistivity  p{T)  is  metallic,  and  p(95  K)  in¬ 
creases  from  110  ^flcm  for  undoped  YBCO  films  to  150 
and  200  /id  cm  in  the  2%  and  3%  Al  doped  films,  respec¬ 
tively.  Scanning-electron  microscopy  (SEM)  pictures 
show  no  evidence  of  misoriented  grains  in  the  surface  to¬ 
pography.  The  most  prominent  defects  we  observe  in  the 
doped  samples  are  holes  about  1000  A  in  diameter  which 
occur  in  groups  of  1-5  clustered  about  1  pm  apart.  The 
clusters  are  typically  20  /xm  apart.  These  holes  may  arise 
from  the  collective  release  of  fluorine  gas  from  the  films 
during  postannealing.  The  similarity  of  our  results  with 
results  on  laser  ablated  films  ^  without  such  defects  indi¬ 
cates  that  they  are  not  of  fundamental  importance  here. 

The  inset  in  Fig.  1  shows  I-V curves  for  the  2  at.  %  Al 
film  at  B  =4  T  parallel  to  c  axis  at  temperatures  evenly 
spaced  by  1  K  from  54.3  to  73.3  K  (all  logarithms  in  this 
work  are  at  base  10).  These  /-Fs  are  typical  of  all  films 
at  all  fields  studied.  We  define  the  glass  temperature, 
T’g=62±l  K,  to  be  the  highest  temperature  at  which  the 
/-K  fits  (to  ±30% )  to  a  single  power  law,  within  our  volt¬ 
age  resolution  (50  nV— 1  mV).  All  /-Fs  above  Tg  exhibit 
systematic  upward  curvature,  thus  restricting  the  uncer¬ 
tainty  of  Tg  to  ±l  K.  Figure  1  summarizes  Tg{B)  [or, 
BgiT)]  for  films  with  0%,  2%,  and  3%  Al.  Note  that  Tg 
for  the  2%  and  3%  Al  films  is  much  lower  than  for  un¬ 
doped  YBCO.  Surprisingly,  doping  at  this  level  weakens 
pinning  and  enhances  flux  wandering.  It  does  not  change 
the  scaling  functions  and  exponents  of  the  transition.  In 
the  inset  of  Fig.  2,  we  demonstrate  how  the  /-Fs  above 
Tg  from  different  samples  coUapse  on  the  same  scaling 
curve.  The  two  fitting  parameters  are  the  low-voltage 
linear  resistivity  p(  T),  and  the  crossover  current  density 
Jq{T)  separating  linear  from  nonlinear  conductivity. 
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FIG.  1.  Transition  temperature  7,  vs  magnetic  field,  for  pure 
YBCO  (circles),  2%  A1  (triangles)  and  3%  AJ  doped  films 
(squares).  The  solid  lines  are  fits  B  with  7*. 

where  p  =  0  at  5  =0.  Inset:  Typical  I-V data  (in  mks  units)  for 
a  2%  AI  film  at  B  =4  T,  from  73.3  to  64.3  K.  The  dashed  line 
indicates  the  7,  at  62 ±  1  K. 


We  extract  the  critical  exponents  z  and  v  of  the  transition 
by  fitting  p  and  Jq  to  the  vortex-glass  relations 
and  Jq  This  provides 

another  way  to  define  7^,  that  leads  to  the  same  value  of 
Tg  (within  2  K)  with  its  operational  definition  from  the 
/-Ts.  Table  I  shows  values  for  the  exponents  v<2  — 1) 
and  2v,  as  well  as  (z  + 1  )/2  from  the  slope  of  the  /-K  at 
Tg  [  which  are  in  excellent  agreement  with 

values  reported  for  YBCO  films^  and  twinned  crystals.^ 
Theoretical  models  (other  than  vortex  glass)  would  pre¬ 
dict  similar  scaling  laws  with  different  expressions  of  the 
exponents  in  terms  of  v  and'z/  In  any  case,  the  fact  that 
scaling  occurs  with  the  same  exponents  and  scaling  func¬ 
tions,  even  as  the  width  of  the  resistive  transition  in¬ 
creases  from  3  to  30  K  with  increasing  doping  and  field, 
supports  the  existence  of  a  phase  transition  at  TgiB)  and 
our  use  of  appropriate  models  to  interpret  our  data. 

Beyond  scaling  arguments,  it  is  possible  to  estimate  the 
magnitude  of  the  vortex  correlation  length  ^c^T)  in  a 


FIG.  2.  Normalized  resistivity  (mks  units)  from  a  3%  Al  film 
vs  r=(  7/7,  - 1 )  for  B  “2,  4,  and  6  T  [with  viz  “  1  )=7.4,  7.2, 
and  7.1  respectively].  The  dashed  lines  indicate  r*  which 
bounds  the  critical  regime  at  high  7.  Inset:  The  scaling  func¬ 
tions  E/Jp  vs  J/Jq  for  a  0%  (lines)  and  a  3%  (circles)  film  col¬ 
lapse  on  the  same  curve. 

model-independent  fashion.  In  particular,  we  arc  in¬ 
terested  in  estimating  temperature  7*  where 

the  scaling  breaks  down,  and  associate  its  value 
to  some  important  mesoscopic  length  scale 
for  vortex  correlations.  We  define  7*  as  the  temperature 
where  p(T)  decreases  5%  below  the  value  predicted  by 
the  scahng  power  law  pa(r-r  To  illustrate, 

Fig.  2  shows  p(  D/( T/Tg  - 1  vs  ( T/Tg  - 1 )  for  the 
3%  Al  doped  film  at  B  =2,  4  and  6  T.  Dividing  out  the 
strongest  T  dependence  of  p  emphasizes  the  break  at  7*. 
The  value  of  T*  and  thus  our  conclusions  arc  not  partic¬ 
ularly  sensitive  to  our  criterion  of  5%,  which  determines 
r*  to  ±  1  K.  Interestingly,  the  critical  scaling  persists 
over  60±5%  of  the  width  of  the  resistive  transition,  even 
though  the  width  varies  from  3  to  30  K.  Here,  we  take 
the  onset  of  the  resistive  transition  to  be  where  p(T)  is 
about  15%  of  the  extrapolated  normal  state  resistivity. 
This  corresponds  to  a  suppression  in  the  mean-field  tran¬ 
sition  temperature  of  ~0.6  K/T,  about  the  same  as  in 
pure  YBCO.  The  significance  of  this  “universal” 
behavior  for  T*  is  unknown. 


TABLE  I.  Critical  scaling  parameters  of  YBa2(CU|_,Al;,  IjOt-j  films.  The  uncertainties  represent 
possible  systematic  errors  in  the  calculations. .  . 


(z  +  l)/2 

V(2-1) 

2v 

t  (±40) 

%  Al 

B  m 

(±0.1) 

(±0.4) 
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We  now  estimate  crossover  current 

density  which  we  define  from  the  scaling  forms,  as 
the  current  density  J  where  E/pJ— 1.5  (inset  of  Fig.  2). 
Since  all  /-Ps  in  the  critical  regime  collapse  on  the  same 
scaling  curve,  other  choices  of  this  ratio  between  1  and  3 
do  not  affect  our  results.  We  use  the  conventional  expres¬ 
sion^*^  Jq  that  does  not  depend  on  the  partic¬ 

ular  model  for  the  phase  transition,  to  determine 
Physically,  this  relation  means  that  nonlinear 
behavior  occurs  when  the  current  density  is  large  enough 
to  produce  translations  of  correlated  vortex  segments  of 
length  |:,|  over  distance  perpendicular  to  the  current,  as 
readily  as  thermal  fluctuations  do.  At  /</^  transla¬ 
tions  of  correlated  segments  are  thermally  induced,  while 
at  />/^,  Lorenz-force-induced  translations  dominate. 
Since  Jq  ^(T the  size  of  vortex  fluctuations 
diverges  at  Tg  as  consistent  with  an  in¬ 

creasing  number  of  vortices  becoming  correlated  and 
moving  coherently  as  F— To  clarify  the  notation,  we 
define  r = ( T /Tg  —  1 )  and  write  = |•(r•  /rV.  Then, 


Figure  3  shows  ^g^-Jq  {T)<f>Q/kT  vs  r/r*  for  a  3% 
A1  film  for  5=2  and  6  T.  The  fact  that  the  two  data  sets 
overlap  means  that  both  the  magnitude  and  scaling  law 
for  T)  are  independent  of  5,  and  hence  the  average  in¬ 
tervortex  spacing.  The  intercept  at  r/r*  =  1  for  this  film 
yields  =  515±40  A,  which  exceeds  somewhat  the  inter¬ 
vortex  spacing  a^=(<f>Q/B)^^^^l90  A  and  320  A  for 
5  =  6  T  and  2  T,  respectively.  Taking  into  account  the 


B(T) 

0  2  4  6 


log  (t/t*) 


FIG.  3.  Plot  of  ^a^=-{J^<l>Q/kT)  of  a  3%  A1  film  for  5  =2  T 
(filled  circles)  and  5  =  6  T  (open  circles)  vs  r/r*.  The  fit  (solid 
line)  corresponds  to  ^^=(515  A)^(r/T*)“^  ®  If  were  pro¬ 
portional  to  the  intervortex  spacing,  then  the  2  T  data  would 
fall  on  the  dashed  line.  Inset:  The  resistivity  p(  T  =  T* )  vs  B 
for  the  undoped  (triangles)  and  3%-Ai-doped  film  (circles).  The 
uncertainty  in  p*  is  indicated  by  the  error  bar  in  one  data  point. 
(All  quantities  in  mks  units.) 


uncertainties  of  our  estimation  in  Tg,T*  and  the  ex¬ 
ponent  v(z-l),  we  find  i.e.,  an  extremely 

weak  dependence  of  on  the  applied  field.  This  is  the 
main  result  of  this  work  and  it  holds  for  all  the  samples 
we  studied,  and  for  the  laser  ablated  films  of  Ref.  1. 
Table  I  records  the  values  of  obtained  for  each  film. 
A1  doping  up  to  3  at.  %  has  a  small  15%  effect  on  the 
magnitude^of  I*.  Evidently  film  morphology  affects 
In  a  250()-A-thick  YBCO  film  with  a  sharp  transition,  we 
found  ^*  =  1240±60  A.  In  an  unusually  resistive  2500- 
A-thick  YBCO  film  with  a  wider  superconducting  transi¬ 
tion  but  the  same  scaling  exponents  we  found 
|'*=732±40  A.  Analyzing  the  measured  Jq{T)  for 
5=2,^  3,  and  4  T  from  Koch  et  aO  on  laser-ablated 
4000-A-thick  YBCO  films,  we  find  |*  =  1100±200  A,  in¬ 
dependent  of  the  magnetic  field. 

What  is  the  significance  of  this  result?  First  we  consid¬ 
er  the  case  that  I**  strongly  depends  on  the  behavior  at 
r*,  which  is  possible  if  T*  marks  the  crossover  be¬ 
tween  two  regimes  with  different  scaling  characteristics. 
A  three-dimensional  XY  critical  scaling  model  proposed^ 
to  describe  the  conductivity  around  is  consistent  with 
our  observation  that  both  boundaries  T*(B)  and  Tg(B) 
scale  with  field  with  the  same  ZY-like  exponent 
[(r,-r*)oc(r,-r^)cc5>/^"',  v'«0.7].  However,  an 
attempt  to  scale  our  p(  T)  data  in  the  vicinity  of  T*  to  the 
expected  form^  cr=ohmic 

conductivity,  coherence  length] 

leads  to  nonuniversal  scaling  functions 
strongly  dependent  on  dopant  concentration.  This  is  true 
even  if  we  consider  an  increase  of  the  magnitude  of 
to  accommodate  the  broadening  of  the  transition  with 
doping. 

We  then  consider  the  magnitude  of  in  terms  of  vor¬ 
tex  correlations.  In  any  model  for  the  phase  transition  at 
Tgy  the  critical  response  must  involve  correlated  motion 
of  multiple  vortices.  Hence,  the  correlation  length 
(or  at  least  should  be  larger  than  the  average  intervor¬ 
tex  spacing  in  the  critical  region.  Models,  in  which  pin¬ 
ning  is  from  quenched  disorder  with  an  average  spacing 
between  pinning  centers  much  smaller  than  the  intervor¬ 
tex  spacing,  lead  to  the  conclusion  that  scaling  should 
break  down  when  the  correlation  length  is  roughly  the  in¬ 
tervortex  spacing  (a  field-dependent  quantity).  In  the 
vortex-glass  theory  in  particular,  the  spatial  extent  of 
phase  coherence  of  the  superconducting  order  parameter 
^G  is  assumed  to  exceed  the  range  of  translational  order 
of  the  flux  lattice.  Thus,  the  critical  scaling  is  expected  to 
break  down  when  equals  the  Larkin-Ovchinnikov 
length,  which  is  at  most  a  few  flux-lattice  spacings.^ 
Again,  should  be  quite  sensitive  to  B  and  T,  in  sharp 
contrast  to  our  observations. 

A  natural  interpretation  of  our  result  is  that  reflects 
the  underlying  pinning  structure  in  the  films.  Quite  gen¬ 
erally,  models  for  critical  scaling  in  samples  where  corre¬ 
lated  disorder  dominates  pinning,  require  the  transverse 
displacement  of  liberated  vortex  segments  to  be  at  least 
the  distance  between  neighboring  macroscopic  defects  in 
the  sample,^  and  thus  lead  to  a  field-independent  value 
for  1“*  at  the  temperature  where  critical  scaling  breaks 
down.  Equally  probable  sites  for  flux  pinning  in  our  films 
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arc  grain  boundaries,  twin  planes  perpendicular  to  the 
substrate,  and  columnar  defects  along  the^  c- 
crystallographic  axis.  A  typical  grain  size  of  500  A  in 
YBCO  films  has  been  measured  by  TEM,  a  value  close 
to  our  Pinning  in  twin  planes  is  also  shown  to  dom¬ 
inate  the  regime  above  the  Tg  in  twinned  YBCO  crys¬ 
tals.  Finally,  assuming  that  only  columnar  defects  such 
as  screw  dislocations  in  crystallites  or  edge  dislocations  in 
grain  boundaries  are  responsible  for  the  pinning^  a  Bose- 
glass^  calculation  leads  to  a  reasonable  d  ^200  A  for  the 
distance  between  adjacent  defects.  Although  our  exper¬ 
iment  cannot  decide  for  the  nature  of  the  extended  de¬ 
fects  responsible  for  the  critical  scaling,  it  suggests  that 
their  structure  is  dictated  by  the  morphology  of  the  sub¬ 
strate  and/or  misfits  in  the  substrate-film  interface  rather 
than  the  dopant  concentration,  since  A1  doping  has  only 
a  15%  effect  on  Furthermore,  the  variation  of 
with  film  thickness  in  undoped  films  suggests  that  defects 
change  somewhat  when  the  film  thickness  exceeds  the 
typical  grain  size.  It  remains  to  be  understood  how  A1 
doping  weakens  the  pinning  properties  of  the  extended 
defects. 

As  a  final  point,  we  discuss  the  striking  result  that  the 
values  of  p*  =p(  T  =  T* )  are  clustered  around  22  /ifl  cm 
for  all  our  samples  and  at  all  fields  (inset  in  Fig.  3). 
Laser-ablated  films  also  show  p^^lS  pilcm  at  fields 
B  =  2—4  T.  *  We  can  show  that  p*  and  set  a  strong 
constraint  for  the  magnitude  of  the  relaxation  time 
t  associated  with  the  critical  slowing  down. 

We  estimate  f  *  by  considering  a  vortex  excitation  of  area 
formed  when  current  density  Jq  passes  through  the 
sample  at  temperature  T*.  The  energy  gained  in  the  sys¬ 
tem  by  the  generation  of  the  excitation  is  Jq  . 

As  this  excitation  relaxes,  it  dissipates  energy  at  a  rate 


f  with  J  and  c  being  the  local  current  density  and 

local  electric  field  caused  by  the  excitation,  integrated 
over  volume.  Assuming  that  the  integral  is  sign^cant 
only  at  the  core  of  the  vortex  of  radius  and  that 
)  is  the  macroscopic  electric  field 
E^P*Jq  increased  by  the  geometric  factor  a^/^aL* 
obtain  f  Letting  f*  be  the  time  for 

the  fluctuation  to  dissipate  its  energy  of  kT,  we  find 

Oisson  et  find  from 
the  complex  impedance  of  a  2500-A-thick  YBCO  film  at 
B  =0.55  T  that  critical  scaling  persists  for  f  >2 X  lO""®  s. 
Using  values  for  and^p*  from  our  2500-A  pure  YBCO 
film,  and  A  at  r*  =  88  K,  we  estimate 

r*»3X  10“*  s,  in  reasonable  agreement  with  the  experi¬ 
ment.  As  B  increases  and  T*  moves  away  from  7^(5), 
inflated  values  of  are  avoided  and  the  system  should 
become  slower  with  a  One  should  note  however, 
that  a  different  picture  for  the  vortex  excitations  could 
lead  to  a  new  combination  of  the  relevant  length  scales 
(including  possibly  an  anisotropic  )  in  the  estimation  of 
t*.  The  issue  may  be  resolved  by  dynamical  measure¬ 
ments  in  an  extended  region  of  magnetic  fields  and  film 
thicknesses.  These  measurements  have  the  potential  to 
estimate  directly  the  relaxation  time  of  the  vortex  excita¬ 
tions,  which  reflects  the  dynamics  of  the  flux  lines  as  they 
cross  and  recombine  in  the  sample. 
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We  present  current-voltage  characteristics  a-V)  in  films  (0^x^0.03),  in  magnetic  fields 

0.5 ^6  T,  which  are  consistent  with  the  existence  of  a  phase  transition  with  critical  exponents  and  scaling 
functions  independent  of  A1  concentration.  We  study  the  role  of  A1  doping  in  flux  pinning  by  analyzing  the  non- 
umversal  characteristics  of  the  transition  such  as  the  width  of  the  critical  regime,  the  characteristic  length  for 
vortex  correlations,  and  the  temperature  and  the  resistivity  that  correspond  to  the  onset  of  the  critical  response. 


1.  INTRODUCTION 

The  nature  of  the  disorder  that  pins  the  flux  lines 
in  high-Tc  superconducting  samples  is  still  a 
controversial  issue.  One  reason  is  that  pinning  of 
the  vortex  system  manifests  itself  through  a  variety 
of  novel  static  and  dynamical  characteristics,  which 
are  quite  different  in  samples  with  varying  degree  of 
disorder  -  as  in  twinned  and  untwinned  crystals,**^  or 
films.^  In  addition,  it  is  shown  theoretically  that 
flux  p  ining  in  extended  or  point-like  defects  may 
lead  to  qualitatively  different  vortex  phases.^  In 
this  work  we  emphasize  the  similarities  in  the 
measured  parameters  in  fV^CO)  films, 

where  we  vary  the  disorder  by  doping  with  Al. 

2.  CRITICAL  SCALING 

It  has  been  demonstrated  that  the  ^c-conductivity 
in  laser  ablated  YBCO  films  can  be  described  by 
appropriate  scaling  laws  that  do  not  depend  on  the 
applied  magnetic  field  B.’  This  fact  supports  the 
idea  for  the  existence  of  a  phase  transition  separating 
a  "liquid*  from  a  *glass"  vortex  phase.  We  test  the 
umversality  of  the  critical  exponents  and  scaling 
functions  at  the  phase  boundary  Tj(B)  through  I-V 
measurements  in  pure  and  Al-dop^  YBCO  films, 
made  by  co-deposition  of  Y,  BaF,  Cu  and  AI,  with 
an  ex-sit u  postanneal  in  oxygen.* 

Typical  I-V  curves  for  a  100  /im  x42  /ira  strip  on 
a  1000  A-thick  film  are  shown  in  the  inset  of  Fig,  1, 
where  we  define  T^  to  be  the  temperature  at  which 
log(Y)  vs.  log(T)  is  closest  to  linear.  In  ref.  7,  we 
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Figure  1.  Scaling  function  forT>T^,  B  =  6T.  Data 
from  a  3  %  (dots)  and  a  0  %  (lines)  film  collapse  on 
the  same  curve.  Inset:  I-Vs  (in  MKS  units)  for  a  2% 
film  at  54.3^T^73.3  K  and  B=4T. 

collected  the  information  from  all  our  films  (two 
pure  YBCO  films,  one  2%  Al  and  one  3%  Al)  in 
the  range  of  fields  0.5  ^B  ^6  T,  and  determined  the 
scaling  exponents  of  the  transition  from  the  I-V  at 
Tj  fV  oc  low-voltage  resistance  [R  oc  (T- 

and  the  crossover  current  separating  the 
linear  from  the  nonlinear  conductivity  oc  (T- 
These  exponents  are  in  excellent 
agreement  with  the  results  from  laser  ablated  YBCO 
films^.  In  Fig,  1,  we  go  one  step  further  to  show 
that  the  scaling  function  describing  the  dc- 
conductivity  above  T^  is  independent  of  Al 
concentration.  All  the  above  argue  in  favor  of  a 


universal  critical  scaling  in  the  vicinity  of  T^. 

3.  TRANSITION  WIDTH 

Doping  YBCO  films  with  up  to  7at%  A1  reduces 
uniformly  the  mean  field  transition  Tg*  by  1.5  K  per 
at%,  and  widens  the  superconducting  transition  at 
B-0.*  In  the  presence  of  a  B  field,  it  drastically 
reduces  T,,  as  shown  in  Fig.  2  (inset).  For 
reference,  our  1%  film  showed  T,=62±0.5  K  at 
B«4T,  to  be  compared  to  T,-76.3  K  for  our 
undoped  film  and  T,»79.5  K  for  a  laser  ablated 
film  in  ref.  3.  This  reduction  can  only  partially  be 
explained  by  the  increase  in  the  penetration  depth 
due  to  doping.  The  energy  gain  per  unit  length 
from  a  vortex  core  pitmed  in  a  defect  t= 
decreases  by  a  factor  of  3  with  2at% 
doping,  given  that  the  zero-temperature  measured 
independently*  increases  to  28(X)  A  relative  to 
1600A  of  pure  YBCO.  However,  T^  decreases  by 
only  20%.  It  is  possible  that  A1  introduces  defecte  in 
the  films  that  tend  to  raise  T,,  like  screw 
dislocations  in  the  crystallites  or  edge  dislocations  in 
the  grain  boundaries.  Such  mechanism  may 

dominate  the  reduction  of  e  at  lower  A1 
concentrations  (x<0.01)  and  consequently  raise  T, 
above  its  value  in  pure  YBCO. 

4.  CORRELATION  LENGTH  AND  TDVIE 

In  Fig.  2  we  illustrate  that  the  resistivity  p(T) 
follows  a  scaling  law  up  to  a  value  of  p‘=22±4  fiQ 
cm,  independent  of  the  dopant  concentration. 
Within  ±3%  error,  p(T)  agrees  with  the  form 

P.(D=P*((T-T^/(r-T^r  (1) 

where  T“  is  the  temperature  at  which 
P(T)“p*.  The  same  value  of  p’  fits  all  our  dau  at 
fields  O.S^B^S  T,  despite  the  fact  that  A1  doping 
doubles  the  normal  state  resistivity  in  samples  with 
3at%  substituted  Cu.  In  addition,  we  find  that  T-T 
is  always  60  ±5  %  of  the  resistive  transition  Tc*-T  * 
and  varies  from  3  to  30  K  for  the  B-fields  and  A1 
concentrations  we  studied.  Doped  samples  may  thus 
be  proved  a  useful  tool  in  case  a  large  critical 
temperature  regime  is  required. 

Eq.  1  sets  a  strong  constraint  for  the  magnitude 
of  the  vortex  correlation  length  {  oc  (T-TJ"  and 


P(T)  (n  m) 

Figure  2.  Deviation  of  the  resistivity  from  power 
law  at  B=4T.  Inset:  Phase  boundary  B=B(Tj). 


time  t  oe  with  i»,  z  the  static  and  dynamical 
critical  exponent  of  the  transition.  In  the  framework 
of  the  vortex  glass  model,  p  oc  (T-Tj)'**  '*  and  thus 
P  «  {/t.  Therefore,  at  T=T,  the  ratio  {/t  is 
independent  (similar  to  p*)  on  the  applied  magnetic 
field,  the  width  of  the  transition  or  the  normal  state 
resistivity  of  the  sample.  The  magnitude  of  {  and  t 
can  be  measured  independently  by  dc~  and  ac-* 
conductivity  measurements  respectively,  and  provide 
quantiUtive  checks  for  the  theoretical  models  that 
describe  the  critical  response  of  the  vortex  system. 
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The  complexity  of  the  properties  of  the  mixed  state  in  the  high-temperature  superconductors 
provides  a  challenge  to  both  fundamental  and  applied  research.  The  central  issue  is  whether 
thermal  fluctuations  can  destroy  the  defining  property  of  superconductors-  a  vanishing  linear 
resistivity.  Early  research  on  the  static  and  dynamical  properties  of  the  mixed  state 
answered  the  question  in  the  affirmative,  suggesting  that  true  critical  currents  could  not  be 
obtained  in  these  materials,  and  dissipation  would  dominate  their  behavior  at  technologically 
useful  temperatures.  However,  as  the  quality  of  the  samples  improves  and  the  sample 
disorder  becomes  better  characterized,  novel  phenomena  are  revealed  indicating  that  there 
is  a  transition  in  the  phase  diagram  from  a  "liquid"  vortex  state  with  finite  linear  resistivity 
to  a  "glass"  state  with  true  critical  currents.  Issues  relative  to  the  nature  of  the  transition 
and  the  role  of  the  disorder  that  pins  the  vortices  in  the  glass  state  are  fiercely  debated  and 
are  not  by  any  means  settled  yet. 

In  our  work  we  focus  on  films,  where  little  has  been  done  towards  recording  the 


1 


2 


statics  and  the  dynamics  in  the  vicinity  of  the  glass  transition  and  in  the  low-temperature 
pinned  vortex  state.  We  demonstrate  that  the  critical  scaling  near  the  transition  does  not 
depend  on  material  properties  of  the  normal  and  superconducting  state,  which  are  altered 
by  doping  the  YBajCujOy^  with  Al.  We  propose  empirical  formulas  that  describe  the 
complex  conductivity  of  the  vortex  phases,  and  we  estimate  the  magnitude  of  the 
fundamental  quantities  of  the  transition,  namely,  the  vortex  correlation  length  and  time.  We 
measure  the  strength  of  pinning  in  the  glass  phase,  correlate  it  to  the  magnetic  penetration 
depth,  and  provide  the  first  evidence  that  the  film  admittance  drops  to  zero  over  a  few 
degrees  Kelvin,  a  strong  indication  for  a  phase  transition  at  the  temperature  where  the 
sample  disorder  becomes  ineffective  to  pin  the  vortex  lattice.  Finally,  we  suggest 
alternative  ways  to  model  the  vortex  excitations  given  the  types  of  disorder  present  in  our 
samples,  and  in  the  light  of  the  theoretical  models  available  at  present. 

This  work  was  supported  by  contract  No.  AFOSR-91-0188. 


PHYSICAL  REVIEW  B 


VOLUME  47,  NUMBER  18 


1  MAY  1993-11 


Absence  of  a  gap  in  the  infrared  reflectance  of  Ni-doped  YBa2Cu307_g  Aims 

Michael  J.  Sumner,  Jin-Tae  Kim,  and  Thomas  R.  Lemberger 
Department  of  Physics,  The  Ohio  State  University,  Columbus,  Ohio  43210 
(Received  2  November  1992) 

We  study  the  effects  of  Ni  doping  on  the  infrared  reflectance  of  YBa2(Cui_;tNijj  )307_6  films  with  up  to 
4%  substitution  of  Cu  by  Ni.  Normal-state  spectra  indicate  that  the  primary  effect  of  Ni  on  the  conduc¬ 
tion  carriers  is  to  increase  their  elastic  scattering  rate.  At  4  at.  %  Ni,  the  elastic-scattering  rate  is  calcu¬ 
lated  to  be  large  enough  to  generate  a  large  absorption  onset  at  the  gap  edge  in  the  superconducting 
state,  if  a  gap  exists.  The  superconducting-state  spectra  show  no  evidence  for  such  a  feature.  Compar¬ 
isons  with  phenomenological  models  for  the  optical  conductivity  show  that  the  conductivity,  and  hence 
the  superconducting  density  of  states,  is  very  gapless. 


L  INTRODUCTION 

A  great  deal  of  effort  has  been  expended  studying  the 
infrared  spectra  of  YBa2Cu307_5,  and  other  oxide  super¬ 
conductors,  as  reviewed  extensively  in  Refs.  1  and  2. 
After  this  effort,  interpretation  of  the  optical  properties 
of  YBa2Cu307«.5  remains  controversial  for  both  the  nor¬ 
mal  and  superconducting  states. 

One  interpretation  is  the  “two-component”  model  in 
which  the  normal-state  optical  conductivity  arises  from  a 
conventional  conduction  band  plus  a  “midinfrared  band” 
involving  bound  electrons  of  uncertain  origin.^  In  this 
model,  the  scattering  rate  1  /r  for  conduction  electrons  is 
roughly  IkT /fi  for  T>Tc,^  and  it  decreases  more  rapidly 
upon  cooling  into  the  superconducting  state.^“^  Hence, 
at  low  temperatures  T <  <T^,  ft/r  is  so  much  smaller 
than  the  superconducting  energy  gap  2A(0)  that  the  ab¬ 
sorption  onset  at  the  gap  frequency  would  be  too  small  to 
observe,  even  for  a  BCS  gap.^  In  this  model,  all  of  the  ab¬ 
sorption  observed  at  10  K  or  lower  is  from  the  midin¬ 
frared  band. 

Another  interpretation  is  the  “one-component”  mod¬ 
el.  The  normal-state  optical  conductivity  arises  from  a 
single  band  of  unconventional  conduction  carriers  which 
are  strongly  coupled  to  an  optically  inactive  boson  of  un¬ 
certain  origin.  The  strong  coupling  generates  an 
inelastic-scattering  rate  that  increases  linearly  with  fre¬ 
quency,  and  thereby  accounts  for  the  anomalously  strong 
absorption  observed  in  the  midinfrared.  The  conduction 
carriers  are  responsible  for  absorption  in  the  supercon¬ 
ducting  state.  Small  features  in  the  reflectance  spectra  at 
150  and  450  cm*”*  are  interpretated  as  gaps  for  excita¬ 
tions  moving  on  the  CuO  chains  and  CUO2  planes,  re¬ 
spectively. 

We  discuss  our  data  in  terms  of  the  two-component 
model,  comprising  Drude  conduction  carriers  plus  a 
midinfrared  band.  Our  conclusions  do  not  depend  on 
this  choice.  Previous  workers*’^  have  applied  this  model 
to  pure  YBa2Cu307_§  and  found  that  the  Drude  conduc¬ 
tivity  corresponds  well  to  the  observed  superconducting 
penetration  depth  nm  and 

normal-state  resistivity  p(r)=m*/neV«0.8  pncm 
( r/K),  where  m  *  is  the  effective  band  mass  and  n  is  the 
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density  of  conduction  carriers.  In  other  words,  when  the 
midinfrared  band  is  subtracted  from  the  experimental 
conductivity  found  by  Kramers-Kronig  analysis,  the 
height  and  width  of  the  remaining  Drude  peak  corre¬ 
spond  to  the  measured  resistivity  and  to  a  scattering  rate 
r-HT)=pm/[fio^HO)l 

We  need  to  make  an  assumption  about  how  the  midin¬ 
frared  band  responds  to  Ni  doping  at  the  4  at.  %  level. 
There  is  as  yet  no  accepted  model  which  would  enable  us 
to  make  physical  estimates.  Hence,  we  appeal  to  the 
spectra  of  Orenstein  et  alJ  on  oxygen-depleted 
YBa2Cu307_5  twinned  crystals.  Interpreting  their  con¬ 
ductivities  at  T <  <T^  as  the  midinfrared  band,  we  see 
that  oxygen  depletion  monotonically  reduces  its  magni¬ 
tude  without  affecting  its  frequency  dependence.  Reduc¬ 
ing  Tc  to  50  K,  by  removing  one-third  of  the  chain  oxy¬ 
gens  to  reach  Og  5,  halves  the  midinfrared  band  conduc¬ 
tivity  for  v<2000  cm‘^  Two  comments  are  in  order. 
First,  analysis  of  our  data  with  a  midinfrared  band  any¬ 
where  in  this  range  leads  to  the  same  conclusions,  pri¬ 
marily  because  the  Drude  conductivity  is  so  much  larger 
than  the  midinfrared  band  below  about  400  cm“'. 
Second,  there  is  no  evidence,  either  from  doping  levels,  or 
Ff,  or  the  Hall  coefficient,  that  4%  Ni  has  an  effect  near¬ 
ly  as  large  as  that  of  removing  33%  of  the  chain  oxygens. 
Therefore,  we  use  for  convenience  the  same  T- 
independent  Lorentz-oscillator  parameterization  of  the 
midinfrared  band  which  Kamaras  et  al.^  used  to  fit  their 
data  on  twinned  YBa2Cu307_g  films.  This  conductivity 
falls  at  the  low  end  of  the  range  found  by  Orenstein 
et  aV 

H.  Ni-DOPED  YBajCujOy-g  FILMS 

Information  about  YBajfCui.^^Ni^ljO^.g  comes  pri¬ 
marily  from  bulk  ceramic  samples.*’’  The  solubility  limit 
is  10  at.  %.  As  in  Zn  doping,  the  crystal  remains  ortho¬ 
rhombic,  suggesting  that  the  CuO  chains  remain  intact 
and  that  a  substantial  fraction  of  the  Ni  goes  into  the 
CUO2  planes.  Smaller  ionic  radius  dopants'®  such  as  Al, 
Co,  and  Fe  make  the  lattice  tetragonal  above  4  at.  %,  and 
apparently  strongly  prefer  chain  sites.  Direct  determina¬ 
tions  are  difficult,  and  site  preferences  for  Ni  remain  con- 
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troversial.  decreases  monotonically,  but  T^ix)  curves 
vary  somewhat  from  group  to  group,  perhaps  because 
different  annealing  procedures  lead  to  different  site 
preferences  for  Ni.  In  our  films,  the  relatively  short 
postannealing  time  of  1  h  at  a  peak  temperature  of 
900  ®C,  compared  with  weeks  at  comparable  temperatures 
for  ceramic  samples,  may  favor  random  siting.  In  the 
end,  Tc(x)  vs  x  for  our  films  falls  among  curves  reported 
for  bulk  ceramic  samples. 

Ni  likely  substitutes  as  Ni^"^  for  Cu^"^,  so  it  should  not 
affect  the  carrier  density  significantly.  The  Hall 
coefficients  Rh^T)  for  3.3  at.  %  and  7  at.  %  Ni  ceramics 
are  identical  to  pure  YBa2Cu307-5  within  the  ±30%  un¬ 
certainty  imposed  by  ceramic  samples.  In  contrast,  R^j 
for  O-depleted  YBa2Cu307-.5  with  r^«66  K  is  about 
four  times  larger  than  for  pure  YBa2Cu307^5.^^  4  at.  % 
Ni  has  a  negligible  effect  on  carrier  density. 

As  discussed  below,  the  dc  resistivity  of  a  patterned  3 
at.  %  Ni  YBa2(Cui^;,Ni^)307«5  film  (Fig.  2)  indicates  a 
residual  resistivity  of  about  250  juHcm,  and  a  corre¬ 
sponding  elastic-scattering  rate  of  about  250  K  X  k  /^.  (k 
is  Boltzmann’s  constant.)  The  resistivity  deduced  from 
our  reflectance  spectra  for  4  at.  %  Ni  are  about  half  of 
this,  which  we  consider  good  agreement. 

The  films  are  made  by  codeposition  of  Y,  BaF2,  Cu, 
and  Ni  onto  (100)  SrTi03  substrates  in  the  presence  of  3 
microtorr  of  O2,  with  a  postanneal  in  moist  oxygen  at 
9(X)®C.*^  The  films  are  nominally  300-nm  thick  and 
about  1  cm  in  diameter.  Thickness  is  determined  from 
measurements  on  identically  prepared  patterned  films,  x 
is  the  nominal  concentration  of  Ni  from  the  measured  Ni 
and  Cu  deposition  rates.  We  have  made  good  films  up  to 
10  at.  %  Ni. 

The  structure  of  each  film  is  characterized  by  x-ray 
diffraction,  scanning  electron  microscopy  (SEM),  and  Ra¬ 
man  spectroscopy.*^  Powder  x-ray  diffraction  indicates 
single-phase  films  highly  oriented  with  the  c  axis  perpen¬ 
dicular  to  the  substrate.  These  results  put  an  upper  limit 
of  about  10%  on  the  volume  of  film  consisting  of 
misoriented  grains  and  second  phases.  SEM  micrographs 
show  surface  structure  on  a  few  percent  of  the  sample 
surface,  examples  of  which  are  in  Ref.  13.  This  structure 
appears  to  be  misoriented  grains  with  dimensions  ranging 
from  0.2  to  5  fim.  These  grains  do  not  show  up  in  Ra¬ 
man  spectra,  even  though  they  should  have  much 
stronger  Raman  signals  than  the  properly  oriented 
grains.  We  conclude  that  misoriented  grains  do  not  ex¬ 
tend  deeply  into  the  film.  Raman  also  is  very  sensitive  to 
some  second  phases  with  large  Raman  signals,  such  as 
BaCu02,  but  reveals  none  of  these  phases  in  the  films. 
The  overwhelming  majority  of  the  volume  of 
our  films  consists  of  highly  oriented  grains  of 

YBa2(Cui-;pNi;^  )307_5. 

Figure  1  shows  the  low-frequency  inductive  response  of 
films  with  up  to  6%  Ni.  The  ordinate  is  the  change  in  the 
inductance  of  a  tightly  wound  10-tum  6-mm  diameter 
pancake  coil  pressed  gently  against  the  sample  film.*^ 
The  maximum  current  density  induced  in  the  film  is 
about  3(X)  A/cm^.  These  data  have  three  important 
features.  First,  the  onset  of  superconductivity  is  very 
sharp,  as  shown  by  the  sharp  “comer”  at  Poorer 


T  (K) 

FIG.  1.  Change  in  inductance  of  a  pancake  coil  pressed 
against  a  film.  The  onset  of  superconductivity  is  sharp,  and 
decreases  smoothly  as  Ni  is  added. 


quality  films  show  a  rounded  onset.  Second,  the  transi¬ 
tion  widths  are  only  a  few  percent  of  ensuring  reason¬ 
able  homogeneity.  Third,  r^(x)  decreases  linearly  with 
increasing  Ni  concentration,  demonstrating  reproducibil¬ 
ity  of  the  fabrication  technique.  The  rate  of  suppression 
of  Tc  with  Ni  doping  is  well  within  the  range  of  values  re¬ 
ported  for  bulk  ceramic  samples.^”'** 

Figure  2  shows  the  resistivities  of  films  made  identical¬ 
ly  to  those  reported  here,  then  patterned  into  strips.*^ 
Our  pure  YBa2Cu307-5  films  have  resistivities  p(100 
K)»1{X)±20  pH  cm;  p{T)  extrapolates  to  zero  at  !r=0, 
and  the  transitions  are  a  couple  of  Kelvins  wide. 
ranges  from  88  to  92  K.  Figure  2  shows  that  3%  Ni  in- 


YBa3(Cu..^Ni^)30,_, 


FIG.  2.  Resistivities  vs  T  of  0%,  3%,  and  7%  Ni-doped 
YBa2(Cui_;^Ni^  )307-.6  films. 
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creases  p(T)  and  causes  an  apparent  residual  resistivity 
p(0)«250  pSlcm,  The  data  at  7%  Ni  continue  this 
trend.  For  comparison,  the  resistivity  of  Zn-doped 
YBa2Cu307^5  twinned  crystals  increases  with  Zn  concen¬ 
tration  at  about  half  the  rate  we  observe  for  Ni  doping.*^ 

We  conclude  that  Ni-doped  YBa2Cu307-5  films  are  of 
sufficient  quality  to  warrant  extensive  optical  measure¬ 
ments.  Although  films  are  twinned  and  less  well  charac¬ 
terized  than  untwinned  crystals  might  be,  and  transmis¬ 
sion  of  some  radiation  through  the  film  into  the  substrate 
precludes  a  Kramers-Kronig  determination  of  the  optical 
conductivity,  still  there  are  several  important  advantages 
of  films  for  optical  measurements.  One  is  that  some 
features  in  the  reflectance  are  traceable  to  the  SrTi03 
substrate,  which  means  that  all  of  the  film  is  probed  by 
the  radiation,  not  just  the  surface  as  for  a  crystal.  Thus, 
the  measured  resistivity  and  inductive  response  of  the 
film  are  those  of  the  carriers  reflecting  the  radiation. 
Another  advantage  is  that,  if  there  is  an  energy  gap,  then 
there  is  an  interference  between  the  reflectances  of  the 
film  and  the  SrTi03  substrate  which  amplifies  the  calcu¬ 
lated  drop  in  reflectance  at  the  gap  frequency.  Both  of 
these  advantages  will  be  exploited  below. 

III.  NORMAL-STATE  REFLECTANCE  SPECTRA 

The  reflectances  of  our  films  are  measured  with  a 
Bomem  Corp.  Fourier  transform  spectrometer  with  a 
spot  diameter  of  7  mm  and  with  an  angle  of  incidence  1 V 
from  normal.  The  high  value  of  the  dielectric  function 
6(6>)  of  YBa2Cu307_5  in  the  infrared  ensures  that  radia¬ 
tion  in  the  film  is  polarized  almost  entirely  in  the  plane  of 
the  film.  The  sample  film  and  a  Au  mirror  are  mounted 
several  centimeters  apart  on  the  cold  finger  of  a  helium- 
flow  cryostat.  The  cold  finger  can  be  raised  and  lowered, 
thereby  placing  either  the  sample  or  the  Au  mirror  in  the 
optical  path  of  the  spectrometer.  Typically,  at  each  tem¬ 
perature  several  alternate  measurements  of  film  and  Au 
mirror  ensure  an  accurate  determination  of  the  relative 
reflectance  of  the  sample.  When  Au  mirrors  are  mounted 
at  both  locations  on  the  cold  finger,  alternate  measure¬ 
ments  produce  reflectances  which  agree  within  1%.  The 
difference  is  reproducible  and  is  taken  into  account  in  our 
analysis. 

Because  some  radiation  penetrates  through  the  films 
into  the  SrTi03  substrates,  a  Kramers-Kronig  analysis 
does  not  produce  accurate  values  for  the  optical  conduc¬ 
tivity  of  the  films.  Therefore,  we  compare  our  reflectance 
data  with  calculated  model  reflectances.  The  complete 
model  for  our  samples  comprises  a  conducting  or  super¬ 
conducting  film  3(X)-nm  thick  on  an  infinitely  thick 
SrTi03  substrate.  This  is  a  good  approximation,  since 
the  actual  substrate  thickness  of  1  mm  is  much  larger 
than  the  longest  wavelength  of  interest  here.  The  spec¬ 
trometer  is  run  at  a  resolution  of  4  cm” ^  to  average  over 
the  resonances  due  to  standing  waves  in  the  substrate. 

The  normal-state  optical  response  of  the  films 
comprises  Drude  carriers  plus  a  broad  midinfrared  band, 
so  that  the  complex  conductivity  is  given  by 

(1) 


The  plasma  frequency  of  the  Drude  carriers  is  taken  to  be 
Vp/c=8200  cm"'  (0^  =  1.55X10**  rad/s;  iuOp  =  l  eV), 
independent  of  temperature  and  Ni  concentration.'  For 
reference,  this  value  of  yields  X(0)=c/du^  =  190  nm, 
and  for  p=  100  p(l  cm,  it  corresponds  to  a  scattering  rate 

1 /T=peo<»p  =2  X 10'*  s"  * =2irc  X 106  cm" ' 

=  mKXk/H.  (2) 

The  midinfrared  band  is  shown  in  Fig.  3.  The  peak 
value  of  the  real  part  is  about  (800  pil  cm)”^  The  lowest 
conduction  carrier  conductivity  at  dc  is  about  (500 
/zfl  cm)”^  for  4%  Ni  at  300  K.  Hence,  both  the  real  and 
imaginary  parts  of  the  Drude  conductivity  dominate  over 
the  midinfrared  band  at  low  frequencies. 

A  four-parameter  model  of  the  measured  conductivi¬ 
ties  and  reflectances  of  SrTi03  at  300  and  85  K  are  shown 
in  Figs.  4  and  5.^^  We  use  the  300-K  data  to  model 
reflectances  at  300  K  and  the  85-K  data  to  model 
reflectances  at  1(X)  and  10  K.  Note  the  strength  and  tem¬ 
perature  dependence  of  the  phonon  at  the  lowest  frequen¬ 
cy;  this  phonon  is  apparent  in  our  spectra  and  is  very  use¬ 
ful  in  determining  the  best-fit  resistivity  of  our  films.  The 
phonons  at  175  and  540  cm”^  are, also  present,  but  weak¬ 
er. 

Let  us  begin  with  pure  YBa2Cu307_5.  Figure  6  shows 
reflectance  spectra  for  50  <v<600  cm”^  of  a  pure 
YBa2Cu307_5  film  at  10,  100,  and  300  K.  The  solid 
curves  are  model  calculations.  Each  data  set  is  a  super¬ 
position  of  spectra  taken  with  different  light  sources, 
beam  splitters,  and  detectors  as  appropriate  for  different 
frequency  ranges.  The  noise  level  varies,  being  especially 
large  near  2(X)  and  550  cm“^  but  is  about  ±0.75%  for 
most  frequencies.  In  general,  the  film  becomes  more 
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FIG.  3.  Real  and  imaginary  parts  of  the  midinfrared  band 
used  to  model  our  films. 


cr(co)  =  (TQA\~-icor)-\-a^i^(co)  . 
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FIG.  4.  Real  part  of  the  conductivity  and  reflectance  of 
SrTi03  at  300  K. 


reflective  as  it  cools,  as  expected  from  a  reduced  resistivi¬ 
ty.  At  10  K,  the  film  is  nearly  perfectly  reflecting  at  the 
lowest  frequencies.  Apparently  the  surface  structure  ob¬ 
served  under  the  SEM  does  not  affect  the  film  reflectance 
here. 

Now  consider  more  detail.  The  300  K  reflectance  does 
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FIG.  5.  Real  part  of  the  conductivity  and  reflectance  of 
SrTiOj  at  85  K. 
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FIG.  6.  Reflectance  of  a  pure  YBa2Cu307_5  film  at  10,  100, 
and  300  K.  Solid  curves  are  model  calculations  with  pdOO  K) 
=  125  fxD,  cm  and  p(300  K)=320  fiO.  cm. 


not  turn  up  toward  unity  at  low  frequency  as  expected  for 
conductors.  Instead,  it  has  a  maximum  at  v«  100  cm“^ 
The  model  calculation,  which  uses  p(300  K)  =  320 
/ill  cm,  reveals  that  this  feature  arises  from  radiation 
passing  through  the  film  and  interacting  with  the  SrTi03 
phonon  at  80  cm“^  At  100  K,  the  feature  is  softened  be¬ 
cause  the  film  is  less  resistive  [the  model  calculation  uses 
p  (100  K)  =  125  /xflcm],  so  less  light  gets  into  the  sub¬ 
strate.  At  10  K,  the  model  calculation  uses  a  gap 
2A=220  cm“'=3.5/:r^  and  p=10  /xflcm.  This  is  the 
clean  limit,  so  the  optical  conductivity  of  the  supercon¬ 
ducting  carriers  is  a  8  function  at  v=0.  The  reasonable 
agreement  justifies  our  choice  of  midinfrared  conductivi¬ 
ty. 

Figure  7  shows  reflectance  spectra  for  a  1%  Ni  film  at 
10,  100,  and  300  K.  The  100  and  300  K  spectra  are 
different  from  those  of  pure  YBa2Cu307_5,  but  not  by 
very  much.  The  model  calculations  use  the  same  resis¬ 
tivities,  p(300  K)  =  320  /ill  cm  and  p(  100  K)  =  1 25  pil  cm. 
The  reflectance  should  be  almost  identical  to  pure 
YBa2Cu307_5  at  such  a  low  doping  level,  so  differences 
indicate  the  “noise”  from  sample-to-sample  variations. 
We  will  see  that  the  gap  feature  expected  from  these  films 
is  much  larger  than  sample-to-sample  variations. 

Figure  8  shows  reflectance  spectra  for  a  2%  Ni  film 
over  a  wide  frequency  range.  The  data  above  1000  cm“^ 
show  that  the  doped  films  behave  the  same  as  undoped 
YBa2Cu307_5  in  this  range,  becoming  steadily  less 
reflective.  Below  600  cm”^  this  film  is  not  quite  as 
reflective  as  pure  YBa2Cu307_5,  and  its  reflectance  does 
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FIG,  7.  Reflectance  of  a  YBa2(Cuo.99Nio.oi)307-6  at  10, 
100,  and  300  K.  Solid  curves  are  model  calculations  with  pdOO 
FC)  =125  ^Hcm  and  p( 300  K)  =  320fiClcm. 
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FIG.  9.  Reflectance  of  a  YBa2(Cuo.96Nio. 04)307 -a 
100,  and  300  K.  Solid  curves  are  model  calculations  with  pdOO 
K)  =  220  fiCl  cm  and  p{ 300  K)  =  520  p£l  cm. 
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FIG.  8.  Reflectance  of  a  YBa2(Cuo.98Nio.o2)307-6  filni  at  10, 
100,  and  300  K.  Solid  curves  are  model  calculations  with  pdOO 
K)  =  200  pVL  cm  and  p(300  K)  =  480  pfl  cm. 


FIG.  10.  Reflectance  ofaYBa2(Cuo.96Nio.o4)307-6  film  at  100 
K  with  model  calculations  for  p=100,  200,  300,  500,  and  700 
pH  cm,  curves  A-E,  respectively. 
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FIG.  11.  RcflcctEncc  of  a  TBE2(Cuo  9^Nio. 04)307— 5  filni  2Lt  300 
K  with  model  calculations  for  p=l00,  200,  300,  500,  and  700 
pn  cm,  curves  A-E,  respectively. 


not  increase  as  much  when  it  cools  from  300  to  100  K. 
Model  calculations,  using  p(300  K)  =  480  /iSlcm  and 
pdOO  K)=  200  pH  cm,  show  that  these  features  indicate 
a  higher  resistivity,  as  expected. 

Figure  9  shows  the  reflectance  of  a  4%  Ni  film.  The 
decreased  reflectances  at  300  and  100  K  relative  to  pure 
YBa2Cu307_5  are  described  well  by  an  increase  in  resis¬ 
tivity  to  p(300  K)  =  520  fia  cm  and  p(100  K)  =  220 
pH  cm.  Let  us  quantify  the  uncertainty  in  the  optical 
determination  of  the  resistivity.  Figures  10  and  1 1  show 
the  100  and  300  K  spectra  for  4%  Ni  with  curves  calcu¬ 
lated  for  several  different  resistivities.  The  most  impor¬ 
tant  feature  in  determining  the  best-fit  resistivity  is  the 
steep  drop  below  100  cm***  because  the  Drude  conduc¬ 
tivity  is  much  larger  than  the  midinfrared  band  in  this 
range.  The  best  fits  are  for  p  =  220±50  pH  cm  and 
520±50  pH  cm,  respectively.  These  values  are  100  and 
200  pH  cm  larger  than  for  pure  YBa2Cu307-.5.  We  con¬ 
clude  that  100±50  pflcm  is  a  reasonable  estimate  of  the 
residual  resistivity  introduced  by  4%  Ni.  This  is  some¬ 
what  less  than  the  250  pH  cm  inferred  from  the  dc  resis¬ 
tivity  of  a  3  at.  %  Ni  film  (Fig.  2),  but  in  reasonably  good 
agreement. 

IV.  SUPERCONDUCTING-STATE 
REFLECTANCE  SPECTRA 

The  spectra  at  10  K  (Figs.  6-9)  show  that  the  films  be¬ 
come  gradually  less  reflective  below  600  cm“^  as  Ni  is 
added.  More  significantly,  no  sharp  features  emerge. 


Since  the  spectra  are  similar,  we  will  concentrate  our 
analysis  on  the  4%  Ni  film  because  it  has  the  largest 
scattering  rate. 

First  of  all,  let  us  calculate  what  would  occur  if  all  of 
the  conduction  carriers  condensed  into  a  conventional  su¬ 
perconducting  state.  We  obtain  the  frequency  depen¬ 
dence  of  the  conductivity  from  the  theory  of  Chang  and 
Scalapino^®  (CS)  for  layered  superconductors  with  weak 
electron-boson  coupling.  In  the  CS  theory,  currents  are 
constrained  to  flow  parallel  to  the  planes  so  the  anoma¬ 
lous  skin  effect  does  not  occur  for  long  mean  free  paths. 
The  CS  model  applies  at  all  temperatures  and  for 
electron-scattering  rates  ranging  from  clean  H/r 
<  <2kTc  to  dirty  ii/r>  >2kTc  limits.  Of  course,  area 
missing  from  the  optical  conductivity  at  finite  frequencies 
is  found  in  a  6  function  at  zero  frequency.  We  examine 
dc  conductivities  of  (7o=l/50,  1/100,  and  1/150  pH  cm, 
corresponding  to  l/r^sSS,  106,  and  160  cm'^  respec¬ 
tively.  The  complete  model  includes  the  midinfrared 
band  and  an  infinitely  thick  SrTi03  substrate. 

Figure  12  shows  the  reflectance  of  the  4%  Ni  film  at 
10  K  together  with  curves  calculated  for  2A=185 
cm”^  =  3.5/:r^  (7’^=76  K)  and  various  resistivities.  Fig¬ 
ure  13  shows  the  same  data  with  curves  calculated  for 
2A=400  cm“^=7.6kr^.  Note  that  the  abrupt  rise  at 
about  70  cm”^  in  the  measured  reflectance  is  associated 
with  a  substrate  phonon  which  behaves  differently  in 
different  substrates.  Also,  the  “break”  at  600  cm“Mn  the 
measured  reflectance  is  associated  with  uncertainty/noise 


FIG.  12,  Reflectance  of  a  YBa2(Cuo.96Nio.o4)307-.6  film  at  10 
K  with  model  calculations  for  2A=185  cm“*^  =  3.5kT.  and 
p=50,  100,  and  200  pflem,  curves  A-C,  respectively.  Abrupt 
changes  in  reflectance  at  70  and  600  cm  “  *  are  artifacts  associat¬ 
ed  with  the  substrate  and  a  change  in  infrared  detectors,  respec¬ 
tively. 
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FIG.  13.  Reflectance  of  a  YBa2(Cuo.96Nio. 04)307-5 
K  with  model  calculations  for  2A=400  cm~^  =  7.6kTc  and 
p=50,  100,  and  200  fiQ,  cm,  curves  A-C,  respectively. 

in  the  pyroelectric  detector  used  above  600  cm“^  In¬ 
creasing  2A  beyond  400  cm"^  moves  the  model  into  the 
clean  limit  in  which  absorption  is  entirely  from  the 
midinfrared  band,  so  there  is  no  gap  feature.  Note  the 
steep  drop  between  200  cm*”^  and  2A(0)  and  the 
minimum  at  2A(0)  in  the  calculated  reflectances.  The 
minimum  results  from  interference  between  the  responses 
of  the  substrate  and  film,  as  discussed  in  Ref.  19.  We  em¬ 
phasize  that  such  strong  features  could  not  be  masked  by 
whatever  is  responsible  for  the  1-2%  discrepancies  be¬ 
tween  the  model  calculations  and  data  in  the  normal 
state.  They  are  absent.  The  conduction  carriers  do  not 
condense  into  a  BCS-like  state  with  a  sharply  defined  en¬ 
ergy  gap. 

The  question  then  becomes,  how  non-BCS  must  the  su¬ 
perconducting  state  be?  We  proceed  by  comparing  physi¬ 
cally  motivated  phenomenological  models  for  the  con¬ 
ductivity  of  the  superconducting  electrons  with  the  data. 
Of  course,  in  each  of  these  models  area  missing  from  the 
optical  conductivity  at  finite  frequencies  resides  in  a  8 
function  at  zero  frequency. 

The  simplest  model  is  motivated  by  the  idea  that  some 
of  the  carriers,  e.g.,  those  on  the  CuO  chains,  remain 
“normal”  at  low  temperatures.  Thus,  the  conductivity  of 
the  conduction  carriers  consists  of  a  fraction  that  go 
superconducting  and  have  a  gap  and  a  fraction  l““/j 
which  remain  normal.  Figure  14  shows  reflectances  cal¬ 
culated  for  /^=0.2,  0.4,  and  0.6  superimposed  on  the  10 
K  spectrum  of  the  4%  Ni  film.  While  the  best  fit  is  for 
/^=0,  keeping  in  mind  the  differences  between  data  and 
model  in  the  normal  state,  we  conclude  that  up  to  40%  of 
the  carriers  could  go  superconducting  with  a  BCS-like 
density  of  states,  if  the  other  60%  remain  normal.  The 
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FIG.  14.  Reflectance  of  a  YBa2(Cuo.96Nio. 04) 307-8  film  at  10 
K  with  two-fluid-model  calculations  for  /,  =0.2,  0.4,  and  0.6; 
2A  =  250  cm~^AJkTc\  and  p=100  pftcm.  The  conductivities 
of  the  normal  and  superconducting  carriers  are  shown  in  Fig. 
15,  labeled  “normal  drude”  and  “CS”,  respectively. 


calculated  curve  for  60%  superconducting  electrons  is 
unacceptable  primarily  because  it  drops  too  steeply  be¬ 
tween  100  and  200  cm”^  and  secondarily  because  its  ab¬ 
solute  value  at  low  frequency  is  higher  than  that  ob¬ 
served.  The  steep  drop  arises  from  the  abrupt  increase  in 
the  real  part  of  the  conductivity  at  the  gap  frequency. 
Any  theory  with  an  abrupt  feature  superimposed  on  a 
smooth  background  conductivity  will  disagree  with  the 
data.  The  best-fit  conductivity  in  this  “two-fluid”  model 
is  very  non-BCS  in  that  at  least  60%  of  the  electrons 
remain  normal.  Note  that  a  factor  of  2.5  reduction  in 
superfluid  density  implies  that  the  penetration  depth  A.(0) 
increases  by  a  factor  of  from  pure  YBa2Cu307_5. 

While  we  are  considering  coexisting  normal  and  super¬ 
conducting  carriers,  we  should  digress  a  moment  to  con¬ 
sider  an  inhomogeneous  film  with  a  putative  normal  layer 
on  top  of  a  superconducting  layer.  Model  calculations  of 
a  three-layer  system  of  normal  film  over  superconducting 
film  over  substrate  yield  essentially  the  same  results  as 
the  above  calculation  because  infrared  radiation 
penetrates  through  both  layers.  Thus,  this  type  of  inho¬ 
mogeneity  could  explain  our  data  only  if  the  supercon¬ 
ducting  layer  were  less  than  40%  of  the  film  thickness. 
Such  a  large  inhomogeneity  is  not  possible  in  our  films. 

Finally,  let  us  examine  phenomenological  models  in 
which  the  optical  conductivity  of  the  superconducting 
carriers  is  smooth,  rather  than  having  an  abrupt  feature 
at  the  gap  edge.  Perhaps  then  we  can  fit  the  data  with  a 
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response  explains  the  effect  of  Ni  on  the  normal-state 
spectra  entirely  in  terms  of  the  dependence  of  the  scatter¬ 
ing  rate  of  the  conduction  electrons  on  temperature  and 
Ni  concentration.  The  relatively  flat  and  featureless 
reflectance  spectra  in  the  superconducting  state  can  be 
understood  in  this  framework  only  if  the  conduction  elec¬ 
tron  response  is  very  gapless  and  therefore  very  non-BCS. 

Experimentally,  we  need  data  on  more  doped  systems 
with  other  dopants  with  different  spin  states  and  site 
preferences.  Measurements  on  doped  crystals  would  be 
very  useful,  for  comparison  in  the  midinfrared  range. 
Also,  we  need  corroborating  evidence  through  direct 
low-frequency  measurements  of  the  penetration  depth,  as 
well  as  measurements  in  the  far  infrared,  where  A.(0)  can 
be  determined  from  the  increase  in  reflectance  toward 


unity  that  must  occur  even  in  very  thin  superconducting 
films. 
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